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Abstract 

We study the statistical mechanics of classical two-dimensional "Coulomb gases" with 
general potential and arbitrary /?, the inverse of the temperature. Such ensembles also 
correspond to random matrix models in some particular cases. The formal limit case 
P = oo corresponds to "weighted Fekete sets" and also falls within our analysis. 

It is known that in such a system points should be asymptotically distributed according 
to a macroscopic "equilibrium measure," and that a large deviations principle holds for 
this, as proven by Ben Arous and Zeitouni [BZj . 

By a suitable splitting of the Hamiltonian, we connect the problem to the "renormal- 
ized energy" W, a Coulombian interaction for points in the plane introduced in [SSI] , 
which is expected to be a good way of measuring the disorder of an infinite configura- 
tion of points in the plane. By so doing, we are able to examine the situation at the 
microscopic scale, and obtain several new results: a next order asymptotic expansion of 
the partition function, estimates on the probability of fluctuation from the equilibrium 
measure at microscale, and a large deviations type result, which states that configurations 
above a certain threshhold of W have exponentially small probability. When f3 — > oo, the 
estimate becomes sharp, showing that the system has to "crystallize" to a minimizer of 
W. In the case of weighted Fekete sets, this corresponds to saying that these sets should 
microscopically look almost everywhere like minimizers of W, which are conjectured to 
be "Abrikosov" triangular lattices. 

keywords: Coulomb gas, one-component plasma, random matrices, Ginibre ensemble, 
Fekete sets, Abrikosov lattice, triangular lattice, renormalized energy, large deviations, crys- 
tallization. 
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1 Introduction 

We are interested in studying the probability law 

(1.1) dF^x 1 ,...,x n ) = ^ J e-^ w ^-' x ^dx 1 ...dx n 

where Zn is the associated partition function, i.e. a normalizing factor such that Pn is a 
probability measure, and where 

n 

(1.2) w n (xi, ...,x n ) = ~y^log \xi - Xj\ + n'Y]V(xi), 
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most often called the Hamiltonian. Here the x^s belong to M? (identified with the complex 
plane C), /3 > is a parameter corresponding to (the inverse of) the temperature and V is a 
potential satisfying some growth and regularity assumptions, which we will detail below. 

The probability law P^ is the Gibbs measure of what is called either a classical "two- 
dimensional Coulomb system" or "Coulomb gas" or "two-dimensional log gas", or "two- 
dimensional one-component plasma", or also "Gaussian /3-ensemble" . It was first pointed 
out by Dyson |Dy| that Coulomb gases are naturally related to random matrices. This is 
somehow due to the fact that e Iog ^*^ log l Xi ~ x i\ is the square of the Vandermonde determinant 
rit<j \ x i ~ x j\ an d thus the l aw PrL m the particular case when V(x) = \x\ 2 and (3 = 2 
corresponds, as shown in [G] (see also |Mej . Chap. 15), to the law of eigenvalues for the Ginibre 
ensemble, which is the set of matrices with independent (complex) Gaussian entries with 
normalized mean and variance. For the general background and references to the literature, 
we refer to the book by Forrester |Foj and references therein. 

The Gibbs measure P^ can also be studied for x%8 belonging to the real line (the one- 
dimensional case). In the context of statistical mechanics (general /3), this corresponds to 
"log gases", and in the context of random matrices {(3 = 1,2,4), to Hermitian or symmetric 
random matrices (whose eigenvalues are always real) . We examine that case in our companion 
paper jSS2j . showing the present study can be extended to handle it. We also point out that 
studying P^ with the Xj's restricted to the unit circle and with (3 = 1, 2, 4 also has a random 
matrix interpretation: it corresponds to the so-called circular ensembles, e.g. in the (3 = 2 
case, eigenvalues of the unitary matrices distributed according to the Haar measure. We also 
plan on examining this case in the future. 

All the above models have been studied quite extensively in the literature, particularly 
from the random matrix point of view (although there one can say the planar or complex 
Ginibre case has attracted less attention than the real or Hermitian case), but also from the 
statistical mechanics point of view, more particularly in the physics literature (particularly 
relevant is e.g. |JLM| ). 

The current research on the random matrix aspect in the complex case focuses on studying 
the more general case of random matrices with entries that are not necessarily Gaussian and 
showing the average behavior is the same as for the Ginibre ensemble, see e.g. [Bat ITVj . This 
is referred to as universality, we are instead limited to exact Vandermonde factors but we 
emphasize that our results are valid for all (3, hence they are not limited to random matrices 
and thus for the proof we cannot rely on any explicit random matrix model. Our results also 
have some universality feature in the sense that they are valid for a large class of potentials 
V. 

The function w n can also be studied for its own sake: it can seen as the interaction energy 
between similarly charged particles confined by the potential V. The case where V(x) is 
quadratic arises for instance as the interaction energy for superconducting vortices in the 
Ginzburg-Landau theory, in the regime where their number is fixed, bounded (see [SS2J, 
Chapter 11). In the case where V is equal to zero on a compact set K and to +oo elsewhere, 
which is not treated here, the minimizers of w n are known as Fekete points or Fekete sets, 
cf. the book of Saff and Totik |ST| for general reference. These are interesting in their own 
right - they arise mainly in polynomial interpolation - and the literature on the question of 
their distribution in various situations is vast. When instead V is a general smooth enough 
function (the situation we treat here), the minimizers of w n are called "weighted Fekete 
points" or weighted Fekete sets, and are also of interest, cf. again [ST]. 
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We will pursue the analysis of these weighted Fekete sets, which can be seen as the formal 
limit /3 — >■ +00 of (jl.ip . in parallel with the analysis of (jl.ip for general f3, and obtain new 
results in both cases. 

In the case of the Ginibre ensemble, i.e. when V(x) = \x\ 2 and j3 = 1, it is known that the 
"spectral measure" v n := ^ Yli=l converges to the uniform measure on the unit disc. More 
precisely P„, seen as a probability on the space of probability measures on C (the spectral 
measures) converges to a Dirac mass at /xo = -lBidx. This is the celebrated "circular law", 
attributed in this case to Ginibre, Mehta, an unpublished paper of Silverstein in 1984, and 
then Girko |Gij . The large deviations from this law was established by Ben Arous and Zeitouni 
|BZj (see Theorem H] below): they showed that a large deviations principle holds with speed 
n~ 2 and rate function 

(1.3) /(„)=/ -k«|»-, 1 *«*<»)+/' M'*(x) 

whose unique minimizer among probabilities is of course the "circle law" distribution ^Ibi dx. 

For the case of a general V and a general /3, we are not aware of a large deviations result 
analogous to |BZj in the literature, however it is not difficult to readapt their proof, with the 
rate function, analogue of (jl,3p . being 

(1.4) /(//) = / - log \x - y\ dn{x) dfi(y) + / V(x)dfi(x). 

Jr 2 xR 2 Jr 2 

Again the spectral measure v n = ^ X^=i converges to the minimizer among probability 
measures of /, called the equilibrium measure, which we will denote ^o- I n the case of weighted 
Fekete sets, the analogue to the circular law has been known to be true for a much longer time: 
it was proved by Fekete, Polya and Szego that - Ya=i converges to the same equilibrium 
measure minimizing / (then also referred to as the electrostatic interaction energy), whose 
description goes back to Gauss, and was carried out with modern mathematical rigor by 
Frostman [FRJ. 

We are interested in examining the "next order" behavior, or that of fluctuations around 
the limiting distribution fiQ. Let us mention that such questions have already been addressed, 
often with the point of view of deriving explicit scaling limits (e.g. [BSH lG]) or laws for certain 
statistics of fluctuations. One can see for example [AHMJ where the authors essentially prove 
that the law of the linear statistics of the fluctuations is a Gaussian with specific variance 
and mean, or also [Rid] for related results. Our approach and results are quite different. 

Recalling / and /io are found through the large deviations at speed re -2 , we look into the 
speed n~ 1 and, while we do not prove a complete large deviations principle at this speed, we 
show there is still a sort of rate function for which a "threshhold phenomenon" holds. This 
analysis is based on an expansion, through a crucial but simple "splitting formula" (which we 
present in Section [LT1 below) of w n (xi, . . . ,x n ), as equal to n 2 I(^o) — § logn plus a term of 
order n, which tends as n — > +00 to the "renormalized energy" W, a Coulombian interaction 
of points in the plane that we introduced in [SSI j and whose definition we will recall below 
in Section 11.21 To be more precise the limit term is the average value of W on the set of 
blow-up limits of the configuration of points x\, . . . ,x n at the scale 1/y/n. It is this average 
that partially plays the role of a rate function at speed n _1 . For a precise statement, see 
Theorem [5j 
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Another way of saying this is in the language of T-convergence (for a definition we refer 
to \Bi\ IDMj . suffice it to know that this is the right notion of convergence to ensure that 
minimizers of w n converge (via their empirical measures) to minimizers of /, i.e. to fj,o): it 
is not very difficult to show (for a short proof, see |SS2j Prop. 11.1) that 2% T-converges as 
n — )■ oo to /, defined in (|1.4p . Here we examine the next order in the 'T-expansion" of w n , i.e. 
we study the T-convergence of ^{w n — n 2 I( / uo)) ) and show that the T-limit is (the average of) 
W. Consequently, after blow-ups at scale y/n, minimizers of w n (i.e. weighted Fekete sets) 
should minimize the (average of the) renormalized energy W. For a precise statement, see 
Theorem [2j 

Before yet giving a precise definition, let us mention that we introduced the renormalized 
energy W in [SSl] for the study of interaction of vortices in the context of the Ginzburg- 
Landau energy of superconductivity (for general reference on the topic, cf. [SS2J ) . Config- 
urations that minimize the Ginzburg-Landau energy with applied magnetic field, exhibit in 
certain regimes "point vortices" that are densely packed (there are n > 1 of them) and are 
expected to arrange themselves in perfect triangular lattices (i.e. with 60° angles), named 
Abrikosov lattices after the physicist who predicted them jA]. The Abrikosov lattices are 
indeed observed in experiments on superconductors^. In |SSlj we made this partly rigorous 
by showing that minimizers of the Ginzburg-Landau energy have vortices that minimize the 
renormalized energy W after blow-up at the scale \/n. The conjecture made in |SSlj . also 
supported by some mathematical evidence (see Section II. 2\i , is then that the minimal value of 
W is achieved by the triangular lattice ; if proven true this would completely justify why vor- 
tices form these patterns. Combining this conjecture with the above conclusion that weighted 
Fekete sets should (after blow-up) minimize W, we thus obtain the conjecture that they also 
should locally form Abrikosov (triangular) lattices. 

We can phrase an analogous question for Coulomb gases or random matrices. Here Pn 
induces a probability measure on the family of blow-ups of (xi, . . . , x n ) around a given origin 
point in E — the parameter of the family — at the scale \/n, a blow-up scale after which 
the resulting points are typically separated by order 1 distances. In the limit n — > oo this 
yields a probability measure on the set of configurations of points in the plane and we may 
ask if, almost surely, the blow-up configurations minimize W. Our results indicate that this 
is not the case, however we are able to prove that there is a treshhold phenomenon, in the 
sense that except with exponentially small probability, the average of W is below a certain 
constant, itself converging to the minimum of W as /3 — >• oo, which indicates crystallisation, 
i.e. if the above conjecture is true, we should see Abrikosov lattices as (3 — > oo. 

To our knowledge, this is the first time Coulomb gases or Fekete sets are connected to 
triangular lattices. 

A corollary of our way of expanding w n is that we obtain a next order estimate of the 
partition function Z%, a result we can already state: 

Theorem 1. Let V satisfy assumptions (jl.lip - (|1.13p below. There exist functions fa, fa 
depending only on V, such that for any [3q > and any j3 > and for n larger than some 
no depending on /3q, we have 

(1.5) n/3/i(/3) < logZf - (-^n 2 I{^) + ^\ogn\ < nf3f 2 (f3), 



1 For photos one can see http://www.fys.uio.no/super/vortex/ 
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with /i,/2 bounded in [/3o,+oo) and such that 

(1.6) lim A (/3) = lim / 2 (/3) = -| 

where a is some constant related to W , and explicited in (jl.36p below. 

This improves on the known results, which only gave the expansion log Zn ~ ^n 2 I(no). 
Let us recall that an exact value for Zn is only known for the Ginibre ensemble case of ft = 2 
and V(a;) = \x\ 2 (see |Mej . Chap. 15): it is Z 2 = n~2 n (" +1 )7r n Ilfc=i an d somewhat lenghty 
computations allow to deduce 

3n 2 n 13 

(1.7) logZ 2 = — — + -logn + n(— 1 + -log 2 + -logTr) + <3(logn) as n -> oo, 

the value of I(no) being | for this potential. On the other hand, no exact formula exists 
for other potentials, nor for quadratic potentials if j3 ^ 2. This is in contrast with the 
one-dimensional situation for which, at least in the case of quadratic V, Zn has an explicit 
expression for every /3, given by the famous Selberg integral formulas (see e.g. [A"GZJ). 

As suggested by the strong analogy between Coulomb gases and interacting vortices in 
the Ginzburg-Landau model, we will draw heavily on methods we introduced in [SSI], such 
as the splitting, blow-up, the use of the ergodic theorem, and particular Ginzburg-Landau 
tools. 

The rest of the introduction is organized as follows: first we give some more notation, give 
the assumptions we need to make on V and state the splitting formula, then we present the 
definition of the renormalized energy and the main results from [SSlj that we will use, and 
finally we state our main results and comment on them. 



1.1 The equilibrium measure and the splitting formula 

We need to introduce some notation, and for this we need to describe the equilibrium measure 
Hq minimizing (|1.4|) among probability measures. 

This description, which is now classical in potential theory (see again [ST] ) says that, 
provided limM-H-oo V(x) — log \x\ = +oo and logV is lower semicontinuous, this equilibrium 
measure exists, is unique, and is characterized by the fact that there exists a constant c such 
that, quasi-every where, 

(1.8) + ^>c and + — = c on Supp(/u ), 
where for any measure /i, denotes the potential generated by fj,, defined by 

(1.9) U'*{x) = - log \x -y\dfi(y) = -2vrA" 1 /i. 

Jm. 2 

Here and in all the paper, we denote by A -1 the operator of convolution by ^~ l°g I " I- h is 
such that A o A -1 = /, where A is the usual Laplacian. We denote the support of by E. 
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Another way to characterize U^ is as the solution of the following obstacle problen^ : It 
is a superharmonic function bounded below by c — V/2 and harmonic outside the so-called 

coincidence set 

(1.10) E = {U^° =c-V/2}. 

This implies in particular that is C ' if V is (see [C]). 

It is now a good time to state the assumptions on V that we assume are satisfied in the 
sequel. 

(1.11) V is C 1 ' 1 and lim V(x) - log |x| = +00, 

\x\—>+oo 

(1.12) In a neighborhood of E, V is C 3 and there exists m, fn > s.t. m < AV^ < m, 



(1.13) V is such that <9£ is C l . 

The assumption on the growth of V is what is needed to apply the results from |ST| and 
to guarantee that (|1.4[) has a minimizer. The other conditions are technical, they are meant 
to ensure that /xo and its support are regular enough, which we will need for example when 
making explicit constructions, and that /j,q never degenerates. Indeed, assumptions (|1.12p - 
(|1.13p . together with Ql.Sp . ()1.9p and the regularity of V, ensure that 

AV(x) 

(1.14) fj,Q = mo(x) dx, where mo(x) = — 1e(x). 

hence for some m, m > we have 

(1.15) m < niQ < m. 

Next, we set C = U 110 + \ — c where c is the constant in (|1.8p . This function satisfies 



(1.16) 



Prom our assumptions on V, in particular assumption (|1.12p . it follows that for some c > 
and x in a neighbourhood of E, 

(1.17) C(x) > c dist(x,£) 2 . 

This follows for instance by using (|1.16p . f 1 1 . 1 2 1) to deduce that AC > r/i outside E, and then 
writing AC in coordinates (r, a) where r is the distance to E, to obtain a lower bound for Crr 
outside but close enough to E (see [KS| . |Frej . |Fri| . [U]). In fact [C] shows that this rate is 
optimal. 

The function C arises in the splitting formula for w n which we now present. As mentioned 
above, expanding the probability density to the next order goes along with blowing-up the 




2 The obstacle problem is a free-boundary problem and a much-studied classical problem in the calculus of 
variations, for general reference see [Fril IKSj . 
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points by a factor y/n. We then denote the blown-up quantities by primes. For example 
x 'i = V^ x i-> m o'( x ') = m o{x), etc .... 

The splitting formula, proven in Section [2] is the observation that, for any x\, . . . ,x n , 



IT 

i=l 



n 1 n 

(1.18) w n (x u ...,x n )= n 2 I{n ) - - logn + -W{-V L H', l K a) + 2n V C(i 

where V -1 := (— d X2 ,d Xl ) and 

/ n N 

(1.19) fl* = -27rA -1 ( - W 



vt=l 



and where, in agreement with formula (|1.27|) below, 



(1.20) W(-V X H', l R a) := lim [ i f \VH'\ 2 + vrnlogr? j . 

r/->0 \^2 JM?\ u v =1 B(x' i ,ri) J 



Letting, for a measure v 



(1.21) F„ M: =(K^ ( " V±ff ' ll2) + 2 "/ C<i '') " ■"•<"■» famES.,** 

[ +oo otherwise, 

the relation (|1.18|1 can be rewritten 

(1.22) w n (xi, ...,i n ) = n 2 I(n ) - ^ logn + nF n ^ J • 

This allows to separate orders as announced since we will see that F n (^2™ =1 is typically 
of order 1. 

We may next cancel out leading order terms and rewrite the probability law (jl.ip as 

(1.23) d&P(xi,...,x n ) = —,e- n 2 F ^i 5 H) dxi . . . dx n 
where 

(1.24) ^:=^ef(" 2/ ^)-f l0 ^). 

As we will see below log Kn is of order n/3, which leads to Theorem [TJ 

1.2 The renormalized energy 

We now define precisely the "renormalized energy" W. In all the paper, Br denotes the ball 
centered at and of radius R, and B(x, R) the ball centered at x and of radius R. 

Definition 1.1. Let m be a nonnegative number. For any continuous function x and o,ny 
vector-field j such that 

(1.25) curl j = 2it{v — to), div j = 
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where v has the form 

(1.26) v = ^^<5p for some discrete set A C M 2 , 

peA 

we let 

(1.27) W(j,X) = lim J / xU| 2 + ^log?/Vx(p) 

^2 7 R 2 VUpeA _ B(p ^ ) ^ 

Definition 1.2. Lei m be a nonnegative number. Let j be a vector field in M. 2 . We say j 
belongs to the admissible class A m □ if (|1.25p . (jl.26p hold and 

(1.28) ^ R ^ is bounded by a constant independent of R > 1. 

For any family of sets {Ur}r>o in K 2 we use the notation x~u R f° r positive cutoff functions 
satisfying, for some constant C independent of R, 

(1.29) \V X u R \ < C, Supp( X u fl ) C V R , X u R (x) = 1 if d(x, U R C ) > 1. 
Definition 1.3. The renormalized energy W is defined, for j G A m , by 

(1.30) W(j) = lim S up Wi ^ X f R \ 

with {xb r }r satisfying (|1.29[) for the family {B R } R>0 . 

In theory, many different j's could correspond to a given v (one can always add the 
gradient of a harmonic function). But as it turns out, they only differ by a constant: 

Lemma 1.4. Let m > and v = ^ peA 5 P , where A C M? is discrete, and assume there exists 
j such that 

(1.31) curlj = lixiy — m), div j = 0, and W(j) < +oo. 

T/ien any oi/ier / satisfying (jl.3ip zs suc/i i/iai j — j' is constant. 
Lf there exists j u such that (|1.3ip holds and such that 

(1.32) Jim / j v = 0, 

R^°°Jb r 

then any other j satisfying (|1.3ip is such that W(j) > W(j v ). 

Proof. Let j, j' be as above, then j — j' can be seen as as complex function of a complex 
variable and from (|1.31|) it is holomorphic. From the finiteness of W(j) and W(j') we deduce 
easily that f B \j — j'\ 2 < CR 2 and using Cauchy's estimate for the coefficients of a power 
series together with a mean value argument it follows that j — j' is constant. 

For the second statement, we deduce from the first statement that j = j u + c for some 
constant vector c 7^ 0, and then 

W(j, XB R ) = W{j u ,XB R ) + C • j iv XB R + ^j- j XB R , 

so that dividing by \B R \, passing to the limit as R — > +00 and in view of (|1.32|) . we find 
W(j) = W(j u ) + \\c\ 2 . □ 



Note that this definition slightly differs from [SSlj : A m here corresponds to Ai-K-m there 
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Note that given u, the above lemma shows that either for all j's satisfying (jl.3ip the limit 
liniR-^oo £ B j„ exists, or it exists for none of them. Both cases may occur. 
The following additional facts and remarks about W are from |SS1] : 

- In the definition (|1.30p . the balls {Br}r can be replaced by other families of (reasonable) 
shapes {Ur}r, this yields a definition of a renormalized energy Wjj, where the letter 
U stands for the family {Ur}. 

- The value of W does not depend on {xb r }r as long as it satisfies (jl.29|) . The corre- 
sponding statement holds for Wu under certain assumptions on the family {Ur}r>o- 

- W is bounded below and admits a minimizer over A\. 

- It is easy to check that if j belongs to A m , m > 0, then j' = -^=j(-/ v /m) belongs to A\ 
and 

(1.33) W(j)=m(w(f)-^logmy 

Consequently if j is a minimizer of W over A m , then j' minimizes W over At- In 
particular 

(1.34) min W = m ( min W log m | . 

V ' An \Ai 2 & J 

- The minimizers and the value of the minimum of Wjj are independent of U. However 
there are examples of admissible j's for which Wjj(j) depends on the family of shapes 
{Ur}#>o which is used. 

- If j G A m then in the neighborhood of p G A we have curlj = 2n(5 p — m), div j = 0, 
thus we have near p the decomposition j(x) = V ± log \x — p\ + f(x) where / is C 1 , and 
it easily follows that the limit (|1.27|) exists if x is compactly supported. It also follows 
that j belongs to L^ oc for any p < 2. 

- Because the number of points is infinite, the interaction over large balls needs to be 
normalized by the volume, as in a thermodynamic limit. Thus W does not feel compact 
perturbations of the configuration of points. Even though the interactions are long- 
range, this is not difficult to justify rigorously. 

- In the case m = 1 and when the set of points A is periodic with respect to some 
lattice "Lu + 7Lv then it can be viewed as a set of n points a±, . . . , a n over the torus 
^(a,v) := M. 2 /CIjU + Tjv) with |T(g^| = n. In this case, the infimum of W(j) among j's 
which satisfy (|1.3ip is achieved by j{ ai } = -V^/i, where h is the periodic solution to 
-Ah = 2Tr(£ i 5 ai -l), and 

(1.35) W(j {ai} ) = ]T G(ai - a,) + vr lim (G(x) + log |x|) 

where G is the Green function of the torus with respect to its volume form, i.e. the 
solution to 

-AG(x) = 2vr ( 5 - t^-^) mT (3j5) . 
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An explicit expression for G can be found via Fourier series and this leads to an explicit 
expression for W of the form Yi^j E{ai — aj) where E is an Eisenstein series (for more 
details see |SS1] and also [BSej ). In this periodic setting, the expression of W is thus 
much simpler than (|1.30|) and reduces to the computation of a sum of explicit pairwise 
interaction. 

- When the set of points A is itself exactly a lattice Zu + TLv then W can be expressed 
explicitly through the Epstein Zeta function of the lattice. Moreover, using results 
from number theory, it is proved in [SSI], Theorem 2, that the unique minimizer of 
W over lattice configurations of fixed volume is the triangular lattice. This supports 
the conjecture that the Abrikosov triangular lattice is a global minimizer of W, with a 
slight abuse of language since W is not a function of the points, but of their associated 
current j {at y 

This last fact allows us to think of W as a way of measuring the disorder and lack of homo- 
geneity of a configuration of points in the plane (this point of view is pursued in [BSej with 
explicit computations for random point processes). Another way to see it is to view W as 
measuring the distance between YpeA^p and the constant m in H . However this distance 
is infinite, both because the domain is infinite and because Dirac masses do not belong to 
H , which is why the actual definition is more involved and requires this "renormalized" 
computation using 77 — > (hence the name, borrowed from [BBH]). 

We may now define the constant a which appears in Theorem [T] and in Theorem [5] below: 

If 1 If 

(1.36) a := — I min W dx = — minW / mo(x) log mo(x) dx, 

7T JE A m Q {x) IT Ai 2 J E 

where we have used (|1.34|) and the fact that, from (|1.14p . J E mo = 1. Note that a only 
depends on V, via the integral term, and on the (so far) unknown constant mirn^ W. 

1.3 Statement of main results 

Our first result identifies the T-hmit of {F n } n , defined in (jl.21[) or (jl .22j) . This in particular 
allows a description of the weighted Fekete sets minimizing w n at the microscopic level. Below 
we abuse notation by writing v n = Ya=i when it should be v n = YaLi ^ «■ For such a 
v, we let v' = Y2=i ^x'. b e measure in blown-up coordinates and j v = —\7 ± H > , where H' 
is defined by (|l,19p — equivalently j u is the solution of div j u = 0, curl j v = 2tv(u' — mo') in 
R 2 which tends to at infinity. (To avoid confusion, we emphasize here that v lives at the 
original scale while j v lives at the blown-up scale.) We also let 



i.e. the push-forward of the normalized Lebesgue measure on E by x i— > (x, j Un (^/nx + •)). It 
is a probability measure onEx L^ c (R 2 ,lR 2 ) (couples of (blow-up centers, blown-up current 
around this center)). 

The limiting object as n — > +oo in the T-limit of w n was v = lim„ u n . In taking the T-limit 
of F n , the limiting object is more complex, it is the limit P of P Un , i.e. a Young measure akin 
to the Young measures on micropatterns introduced in |AM| . Note that in the rest of the 
paper, the probability P has nothing to do with P„, and depends on the realizations of the 
configurations of points. 



(1.37) 
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We will here and below use the notation 
(1.38) D{x',R)=v n (B\x,^y\ -niM>\B\x,-j= 

where x' = \Jrix as usual. 

Theorem 2 (Microscopic behavior of weighted Fekete sets). Let the 'potential V satisfy as- 
sumptions (fTTTD - (fTT3D . Fix 1 < p < 2 and let X = E x Lf oc (IR 2 ,lR 2 ). 

A. Lower bound. Let v n = Yli=i^i be a sequence such that F n (u n ) < C. Then P Un 
defined by (|1.37p is a probability measure on X and 

1. Any subsequence of {Pv„} n has a convergent subsequence converging to a probability 
measure on X as n — > oo. We denote by P such a limit. 

2. The first marginal of P is the normalized Lebesgue measure on E. P is invariant by 
(x,j) i y (x,j(X(x) + •)), for any X(x) of class C 1 from E to M? (we will say T^ x y 
invariant ). 

3. For P almost every (x,j) we have j G A mQ ( x y 
4- Defining a as in (|1.36p . it holds that 



(1.39) liminf \F n (u n ) - 2 J Qdv^j > M J W(j) dP{ 



x,j) > a. 

B. Upper bound construction. Conversely, assume P is a T\i x \ -invariant probability 
measure on X whose first marginal is -^dx\E and such that for P -almost every (x,j) we have 
j G A mo ( x y Then there exists a sequence {u n = Ya=i ^iln °f empirical measures on E and 
a sequence {j n }n in -^f Q c(^ 2 '^ 2 ) suc ^ that curlj n = 2i:{y , n — mo') and such that defining P n 
as in (|1.37p , with j n replacing j Un , we have P n — > P as n — >• oo and 



(1.40) limsupF n K) < l^i I W K (j)dP(x,j), 

n— >oo 7T J 



where Wk is the renormalized energy relative to the family of squares {Kr = [— R, R] 2 }ji. 

C. Consequences for minimizers. If (x%, . . . ,x n ) minimizes w n for every n and v n = 
Sr=i ' then ^e limit P of P Un as defined in (| 1.37ft satisfies the following. 

1. For P- almost every (x,j), the current j minimizes W over A mo t x y 

2. We have 

\E\ I n 
lim F n (u n ) = l —± / W(j) dP{x,j) = a, lim V dist 2 (^, £) = 0. 

J i=l 

3. There exists C > such that for every x' G M 2 , R > 1 and using the notation (|1.38p we 
have 

(1.41) D(x', Rf min (l, < Cn. 
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Note that part B of the theorem is only a partial converse to part A because the constructed 
j n need not be divergence free, hence in general j n ^ j Un . Also the lower bound (jl.39p uses W 
and the upper bound (|1.40p uses Wk- In fact (|1.39p is true with W\j for any family {\Jr}r 
satisfying (|6,2|) below, and in particular we could replace W by Wk there. We recall though 
that both have the same minimizers and minimum value, as proved in [SSlJ. 

This theorem is the analogue of the main result of |SSlj but for w n rather than the 
Ginzburg-Landau energy. It is technically simpler to prove, except for the possibility of a 
nonconstant weight mo(x) which was absent from [SSlj. It can be stated as the fact that 
— J WdP, which can be understood as the average of W with respect to all possible blow-up 
centers in E (chosen uniformly at random), is the T-limit of w n at next order. Its minimum 
over all admissible probabilities is a. 

The estimate (|1.41|) gives a control on the "discrepancy" D (between the effective number 
of points and the expected one) at the scale R/y/n. Note that in a recent paper [AOC| . 
the authors also study the fine behavior of weighted Fekete sets. Using completely different 
methods, based on Beurling-Landau densities and techniques going back to |La| . they are able 
to show the very strong result that 

V V D ( X n,R) _ 

lim sup lim sup — ^ = U, 

i?— >oo ri— >co R 

as long as dist(x n , <9£") > log 2 n. This shows that the density of points follows /Xq at the 
microscopic scale 1 / \/n and thus the configurations are very rigid. This still leaves however 
some uncertainty about the patterns they should follow. On the contrary, our result is less 
precise about D{x,R) since we only recover the optimal estimate when R grows faster than 
n 1 / 4 , but it connects the pattern formed by the points to the Abrikosov triangular lattice via 
the minimization of W. 

We now turn to Coulomb gases. It is straightforward from the form (I1.23P and the estimate 
(provided by Theorem [T]) log Kn = 0(n/3) where Kn is defined in (|1.24|) . to deduce that 
F n < C except on a set of small probability. This fact allows to derive various consequences, 
the first being estimates on the probability of certain rare events. 

Theorem 3. Let V satisfy assumptions f)l. 11 j) fjl. 13[) . 

There exists a universal constant Rq > and c, C > depending only on V such that: For 
any (3q > 0, any n large enough depending on /3q, and any f5 > (3q, for any x\, . . . ,x n E M 2 , 
any R > Rq, any x' = \/nxo G M 2 and any rj > 0, letting u n = Y^i=x&xi> we have the 
following: 

(1.42) logpP(\D(x' ,R)\ >riR 2 ) < -cj3 min(r/ 2 ,if)R A + C '(3(R 2 + n) + Cn. 

(1.43) logP^ (J Qdv n >r^j< ~n/3r] + Cn(fi + 1). 
Moreover, for any smooth bounded U' = \fnU C M 2 , 



(1.44) tosK(f 



jD(x '' R)2 min ( 1, \ D ^ R ^ ) dx'> V )< nfii-cr, + C\U\ + C) + Cn. 



R 2 \ ' R 2 
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Finally, if q £ [1,2) there exists c, C > depending on V and q such that V77 > 1, R > 0, 
(1.45) logP^l + ^ 2 V-«wllw-M(B fl/ ^)>W^ < -cn/V + Cn(/3 + 1), 

where W~ l,q {Q.) is the dual of the Sobolev space Wq' p (Q) with | + | = 1; in particular W^ 1,1 
is the dual of Lipschitz functions. 

These estimates can roughly be read in the following way: as soon as rj is large enough, 
the events in parentheses have probability decaying like e~ cn . More precisely, we bound 
the probability that a ball contains too many or too few points compared to the expected 
number n/j,o(B), but whereas the circular law does it for a macroscopic ball, i.e. for R 
comparable to y/n, the estimate (|1.42p is effective at microscopic scales, of the order of n 1 / 4 . 
This is sometimes called in this context "undercrowding" or "overcrowding" of points, see 
|JLM[ INSV1 IK] . In view of similar results in [JLM| and the result of |AOCj , we could expect 
this to hold as soon as R 3> 1, but this seems out of reach by our method. This can also 
be compared with analogous estimates without error terms proven in the case of Hermitian 
matrices, cf. [ESY| . These results, in the Hermitian case, are proven in the general setting of 
Wigner matrices, i.e. Hermitian matrices with random i.i.d. entries, which do not need to be 
Gaussian. They only concern some fixed /3 however. 

The estimate (|1.44p gives a global version of this result: it expresses a control on the 
average microscopic "discrepancy" D. This control is in L 2 for large values of the discrepancy, 
and in L 3 for small values. The estimate (|1 .43|) allows, in view of (|1.17p . to control (again 
with some threshhold to be beaten) the probability that some points may be far from the set 
E. Note that since v n is a non-normalized empirical measure, (jl.43p ensures for example that 
the probability that a single point lies at a distance 77 from E is exponentially small as soon 
as r] is larger than some constant. All these estimates rely on controlling D and F n by W . 

Finally, (|1.45p tells us that fluctuations around the law nfJLo can be globally controlled (take 
for example R = \fn) by 0{\fn) (except with exponentially small probability). We believe 
this estimate to be optimal. Its proof uses in a crucial manner the result of [STJ, which 
controls, via Lorentz spaces, the difference v n — nfiQ in terms of the renormalized energy W. 

Our last result mostly expresses Theorem in a "moderate" deviations language. Before 
stating it, let us recall for comparison the result of [BZ]: 

Theorem 4 (Ben Arous - Zeitouni). Let (3 = 2 and V(x) = \x\ 2 . We denote by P« the image 
of the law (jl.ip by the map (x\, . . . ,x n ) h-» v n , where v n = ^ Y17=l <W Then for any subset 
A of the set of probability measures on M 2 (endowed with the topology of weak convergence), 
we have 

- inf T(p) < liminf ^logP^(yl) < limsup -=■ logP^(A) < - inf T(p), 
where 1 = 1 — min I. 

Before stating our theorem, we introduce the following notation, which allows to embed 
C n into the set of probabilities onI = £x Lf oc (IR 2 , R 2 ): For any n and x = {x\, . . . , x n ) £ C" 
we let i(x) = P Un , where v n = 'Y^l=\ and P Un is as in (|1.37p . so that i(x) is an element of 
7{X), the set of probability measures onI = Bx Lf oc (R 2 , 1R 2 ). 
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Theorem 5. For any n > let A n C C n . Denote 



(1.46) Aoo = p| (J i(A m ). 

n>0 m>n 

Then for any n > there is C v > depending on V and r\ only such that a being as in (|1,36D , 

(1A7^ v logP^AJ P(\E\ . f f wr ,,j,, C v 

(1.47) hmsup < ml / W {j)dP{x, j) — a — n - 1 

n->oo n 2 V 7r pgA^ J p 

Conversely, let A C ^{X) be a set of TW) -invariant probability measures on X and let A 
be the interior of A for the topology of weak convergence. Then for any n > 0, there exists a 
sequence of subsets A n C E n such that 

P f\ E \ ^ fw.,,U D ^ * „. , „ , C V\ ^^.^^n(An) 



fl.48) --(^inf fwK(j)dP(x,j)-a + ri + ^- ) <liminf 

2 \ TT Pe AJ P ) n ^°° 



n 



and such that for any sequence {v n = Ya=i $xi}n such that (x\, . . . ,x n ) G A n for every n 
there exists a sequence of currents j n £ L^ oc (M 2 ,R 2 ) such that curlj n = 1^(y' n — mo') and 
such that — defining P n as in (|1.37p with j n replacing j Un — we have 

(1.49) limine i. 



Note that if P n was P Un , then (|1.49p would be equivalent to saying that C] n L) m>n i(A m ) C A. 
The difference betwen P Vn and P n is that the latter is generated by a current j n which is not 
necessarily divergence free. 

Compared to TheoremS]this result can be seen as a next order (n instead of n 2 ) deviations 
result, where the average of W over blow-up centers plays the role of a rate function, with 
a margin which becomes small as f3 — > oo. While Theorem [3] said that empirical measures 
at macroscopic scale converge to /io, except for a set of exponentially decaying probability, 
Theorem [5] says that within the empirical measures which do converge to fio, the ones with 
large average of W (computed after blow-up) also have exponentially decaying probability, 

2 

but at the slower rate e~ n instead of e~ n . More precisely, there is a threshhold C//3 for some 
C > 0, such that configurations satisfying 

\e\ r c 

l —t / WdP >a + - 
vr J P 

have exponentially small probability, where we recall a is also the minimum possible value 
of M/ WdP. Since we believe that W measures the disorder of a (limit) configuration of 
(blown up) points in the plane, this means that most configurations have a certain order. The 
threshhold, or gap, C / j3 tends to as f3 tends to oo, hence in this limit, configurations have 
to be closer and closer to the minimum of the average of W, or have more and more order. 

Modulo the conjecture that the minimum of W is achieved by the perfect "Abrikosov" 
triangular lattice, this constitutes a crystallisation result. Note that to solve this conjecture, 
it would suffice to evaluate a, which in view of Theorem [1] is equivalent to being able to 
compute the asymptotics of Z% as /3 — > oo. More generally, the following open questions 
naturally arise in view of our results, and are closely related to one another: 

- Prove that min_4 1 W is achieved by the triangular lattice. 



14 



- Find whether a large deviations statement is true at speed n, and if it is, find the rate 
function. 



- Describe the limit of the probability measures P« on T(X) defined as the images of P« 
by the embedding i. 

In [SS2] we will see that all the results we have obtained here are also true in the case 
of points on the real line, i.e. for ID log gases or Hermitian random matrices. There the 
minimization of W is solved (the minimum is the perfect lattice Z) and the crystallisation 
result is complete. 

The rest of the paper is organized as follows: Section 2 contains the proof of the "splitting 
formula". In Section 3, we present the "spreading result" from |SSlj and some first corol- 
laries. In Section 4, we present an explicit construction which yields the lower bound on 
Zn, whose proof is postponed to Section 7. In Section 5, we show how W controls the over- 
crowding/under crowding of points, and prove Theorem[3l In Section 6 we present the ergodic 
averaging approach (the abstract result) and apply it to conclude the proofs of Theorem [21 [5] 
and [Q 

Acknowledgements: We are grateful to Alexei Borodin, as well as Gerard Ben Arous, Amir 
Dembo, Percy Deift, and Alice Guionnet for their interest and helpful discussions. E. S. was 
supported by the Institut Universitaire de France and S.S. by a EURYI award. 

2 Proof of the splitting formula 

The connection between w n and W originates in the following computation 
Lemma 2.1. For any x±, . . . , x n and letting v = 5 Xi the following holds 



1 n 
(2.1) F n (u) = — W(-V ± F',l R2 ) + 2VC(x i ) 

*nir — * 



mr 

1=1 



= - (w n (xi, ...,x n )- n 2 I(fi ) + ^ logra 
n \ 2 

where F n is defined in (|1 .21[) , W is defined in (jl.27p , and H' is defined in (|1.19p . 

Proof. Let v = 2~2i=i <W First we note that since v and n//o have same mass and compact 
support we have H(x) = 0(l/\x\) and VH(x) = 0(l/|x| 2 ) as \x\ — > +oo. 
We prove that, denoting A the diagonal in M 2 x R 2 , 

(2.2) / -log\x-y\d(u-n l xo)(x)d(u-n^ )(y) = -W(-V ± H,l R 2). 

J(R 2 xR 2 )\A 

First, using Green's formula, 



B R \UiB(xi,ri) JdB R du ^r 1 J dB(xi,n) du J B R \UiB(xi,ri) 

Let Hi(x) := H{x) + log|x — x%\. We have Hi = — log*(fj — n/^o), with V{ = v — S Xi , and 
near Xi, Hi is C . Therefore, using p.!9p and the boundedness of fiQ in we have that, 
as rj — y 

dH 

H— = -27Tlogn + 27rHi(xi) + o{l), 



L 



8B{xi,r)) 
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while the integral on <9-Br tends to as R — > +00 from the decay properties of H. We thus 
obtain, as r\ — > and R — > +00, 

/ \VH\ 2 = -2vrnlog?7 + 27r - 2vrn / Hd/i + o(l), 

JB R \i>iB(xi,ri) i J 

and therefore, by definition of W, 

(2.3) W(- V^H, l R 2)=TrJ2 H i( x i)-™J H d Mo- 

i 

Second we note that given x £ M?, either x ^ {xi} and 

/ -log|a;-y|d(^-n^o)(y) = H(x), 

JR 2 \{x} 

or x = x,; and 

log \x - y\ d{v - nfi )(y) = Hi(xi). 



I 

JR 2 \{x} 



It follows that 



J - log \x - y\ d{y - nfj, )(x) d(v - nfj, )(y) = ^ Hi(xi) - n j H(x) dfj, Q (x), 



which together with (J2T3J) proves (|2.2p . 

On the other hand, we may rewrite w n as 

w n (x\, . . . , x n ) = J — log \x — y\ dv(x) dv(y) + n J V{x)dv(x) 
and, splitting u as n/^o + v — n^iQ and using the fact that /Uo(A) = 0, we obtain 
w(xi, . . . , x n ) = n 2 I(po) + 2n J (x) d(v — nfio){x) + n J V(x) d(v — nfj,o)(x) 

+ / - log \x - y\d(v - n no) (x)d(v - n[i )(y). 

Since + \ = c + £ and since v and ra/io have same mass n, we have 



2n I U^{x)d{v -n^){x)+n J V {x) d{v - n^){x) = 2n J ( d{v - fj, ) = 2n I C,dv, 
using the fact that C = on the support of ^o- Therefore, in view of (|2.2p we have found 
(2.4) w( Xl , ...,x n )= n 2 I(p ) + 2nJ(dv + ^W{-V L H, l K a). 

But, changing variables, 



2 iR 2 \Uf =1 B(x ijT ?) 2 JR 2 \U'^B(x' i ,y/E v ) 

and by adding 7rnlog?7 on both sides and letting 77 — ^ we deduce that W(— V^H, 1 R 2) = 
W{-V L H', 1 R2 ) - § nlogra. Together with (|23J) this proves ([13]). □ 
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3 A first lower bound on F n and upper bound on 

The crucial fact that we now wish to exploit is that, even though W(J, x) does not have 
a sign, there are good lower bounds for F n . This follows from the analysis of [SS1| . more 
specifically from the following "mass spreading result", adapted from |SSlj . Proposition 4.9 
and Remark 4.10 (with slightly different notation), which itself is based on the so-called 
"ball construction method", a crucial tool in the analysis of Ginzburg-Landau equations. 
This result, that we will use here as a black box, says that even though the energy density 
associated to W(j, x) is not positive (or even bounded below), it can be replaced by an 
energy-density g which is uniformly bounded below, at the expense of a negligible error. 
For any set Q, n denotes its 1-tubular neighborhood, i.e. {x £ M 2 , dist(x, n) < 1}. 

Proposition 3.1. Assume $7 C R 2 is open and {v,j) are such that v = ^Y^p^A^p f or some 
finite subset A of n and curlj = 2n{i> — m), div j = in n, where m E L°°(n). Then there 
exists a measure g supported on n and such that 

- Given any constant < p, there exists a family of disjoint closed balls B p covering 
Supp{u), with the sum of the radii of the balls in B p intersecting with any ball of radius 
1 bounded by p, and such that 

(3-1) g> -C(\\m\\ L ~+l) + ^\j\ 2 l n \ Bp inn, 

where C is a universal constant. 



(3.2) g = i|j| 2 outside U p& aB{p,C) 

where C is universal. 

For any function x compactly supported in n we have 



(3.3) 



< CN(\ogN+ |MU«)||Vx||c 



W(j,x) - J Xdg 

where N = G A : B(p, C) f] Supp(Vx) ^ f or some universal C . 
For any U C n, 

(3.4) #(A nU)<c(l + \\mf LOB \U\ + g(U)) . 



Note that the result in [SSI] is not stated for any p but a careful inspection of the proof 
there allows to show that it can be readapted to make p arbitrarily small. 

Definition 3.2. Assume v = Y17=i^ x i- Letting u' = ^27=l^x' be the measure in blown-up 
coordinates and j u = — V H' , where H' is defined by (|1.19p . we denote by g v the result of 
applying the previous proposition to (v',j u ) in IR 2 - 

Even though we will not use the following result in the sequel, we state it to show how we 
can quickly derive a first upper bound on Zn from what precedes. 
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Proposition 3.3. We have 

(3.5) log Kl < Cn(3 + n(log \E\ + o(l)) 
and 

(3.6) logZf < -^n 2 I( M ) + ^logn + Cn/3 + n(log| J E| +o(l)) 

where o(l) — )■ as n — > oo uniformly with respect to f3 > /?o, /or any /3o > 0, ana" C depends 
only on V. 

The proof uses two lemmas. 

Lemma 3.4. For any v = Y11=i we have 

(3.7) F n {u) = — [ dg u + 2 [ (dv, 

ravr J R 2 J 

where F n is as in (|1.2ip . 

Proof. This follows from (|3.3[) applied to xb B j where is as in (|1.29p . If R is large enough 
then #{p G Supp(i^) : B(p, C) n Supp(Vx) ^ 0} = and therefore (|3.3|) reads 

w^o^jxbh) = y dpi/. 

Letting R — > +oo yields W(j u , 1 R 2) = J and the result, in view of (|1.2ip . □ 

Lemma 3.5. Letting u n stand for Y17=i ^ we have, for any constant A > and uniformly 
w.r.t. f3 greater than any arbitrary positive constant (5q, 

(3.8) lim ( [ e-^ n ^ dv -d Xl ...dx^] n =\E\. 

n ^°° \Jc n J 

Proof By separation of variables, we have 

e -A/3nE? =1 Cfe) dx x ...dx n =( [ e -^<(-) dx 



On the other hand, we have £ > and = 0} = E by (|1.16p . hence we have e —^ n Cm \ E 
pointwise, as fin — > oo. By Lebesgue's dominated convergence theorem, it follows that (|3.8p 
holds uniformly w.r.t. /? > /3o, for any /?o > 0. □ 

Proof of Proposition [33 Let again f n stand for J^iLi ■ From (|3.2|) we have g Vn > outside 
UiB(xi,C) and from ()3.ip we have g Un > —C (depending only on H/ioIlL 00 hence on V) in 
UiB(xi,C). Inserting into (|3.7|) we deduce that 

F n (v n ) >~C + 2j (dv n , 

where C depends only on V. Inserting into (jl.23p and integrating over C n , we find 

1 < \e Cn P [ e -™/3K^ d Xl . ..dx n . 



Inserting (|3.8|) and taking logarithms, it follows that 

logA'f < Cn/3 + n(log \E\ +o(l)). 
The relation (f3T6|) follows using ([04]) . □ 
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4 A construction and a lower bound for 

In this section, we construct a set of explicit configurations whose W is not too large, and 
show that their probability is not too small, which will lead to a lower bound on Zn- This 
is the longest part of our proof. The method is borrowed from |SSlj but requires various 
adjustments that we shall detail in Section \7\ We will need (|1.13j) in order to simplify the 
construction and estimates near the boundary. 

We recall the following result ( [SS1| . Proposition 4.2). In this proposition, the notation 
Wk is used for the renormalized energy relative to the family of squares Kr = [-R, R] 2 , i.e. 

defined by Wx(j) = lim sup^^ — where Xk r satisfies (|1.29p relative to the family of 

centered squares. 

The following proves Theorem [21 part B and contains a bit more information useful for 
proving Theorem [5] 

Proposition 4.1. Let P be a T\r x \-invariant probability measure on X = E x L^ oc (M 2 , M?) 
with first marginal dx\E/\E\ and such that for P almost every (x,j) we have j G A mo ( x y Then, 
for any i] > 0, there exists 5 > and for any n a subset A n C C n such that \A n \ > n\(n5 2 /n) n 
and for every sequence {u n = 6 yi + ■ ■ ■ + 5 yn } n with (y%,..., y n ) £ A n the following holds. 

i) We have the upper bound 

(4.1) limsup- (w n (yi,. . . ,y n ) - n 2 I(fi ) + 77 log n) < — [ W K (j) dP(x,j) + rj. 

n->oo n \ 2 J 7T J 

ii) There exists {j n }n m L^ c (M 2 , R 2 ) such that curlj n = 2TT{v' n — mo') and such that the 
image P n of dx\ E /\E\ by the map x i-> {x, j n (yjnx + •)) is such that 

(4.2) limsupdist(P n ,P) < rj, 

n— >oo 

where dist is a distance which metrizes the topology of weak convergence on y{X). 

Applying the above proposition with 7] = 1/k we get a subset A n ^ in which we choose 
any n-tuple (jJi k)i<i<n- This yields in turn a family {P n ^k} of probability measures on X. A 
standard diagonal extraction argument then yields 

Corollary 4.2 (Theorem [21 Part B). Under the same assumptions as Proposition ^. 1[ there 
exists a sequence {v n = S Xl + ■ ■ ■ + d~ Xn } n and a sequence {j n }n m -^f oc 0^ 2 > suc h that 
curl j n = 2-K(v' n — mo') and 

(4.3) limsup - (w n (xi, ... ,x n ) - n 2 I(^ ) + 77 logn) < — [ W K {j) dP(x, j). 

Moreover, if we assume in addition that P is Tu x \ -invariant then, denoting P n the image of 
dx\E/\E\ by the map x h-» (x,j n (yjnx + ■)), we have P n — ^ P as n — ^ +00. 

Another consequence of Proposition 14.11 is. recalling (jl.36p and ()1.24D : 

Corollary 4.3 (Lower bound part of Theorem [1]). For any r) > there exists C v > such 
that for any f3 > we have 

(4.4) hminf— > --(a + r/) - C v . 

n 2 
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Proof of the corollary. We use Proposition 14.11 Choose jo £ to be a minimizer for 
min^j W, and let P be the image of the normalized Lebesgue measure on E by the map 

x i-> (x,a mo(x) j ), where 

(4.5) o- m j(y) := xjm j \\fmy) . 

Then by construction P-almost every (x,j) satisfies j E ■A. mo { x ) an d the first marginal of P is 
dx\ E /\E\. 

Given n > 0, applying Proposition 14.11 and using the notation there we have \A n \ > 
n\{5 2 /n) n and from (|1.23p we have 

(4.6) 1 > f * e -4^K) dyi . . . dynj 

Ja„ K n 

where v n = Yn=i <W F rom f|l-22j) and (|4.ip . when (yx, • • • , 2/n) £ we have 



f„(i/„) < v + — [ w K {j) dP(x,j) = n + - f w K (c 

k J vr J E 



'«() 



From (jl.33p . which is true for Wk as well as W, and since f E mo = 1, we obtain 

- / W K {<y mo ( x )jo) dx = -W K {jo) ~\ [ m (x) log m (x)dx = a, 
n Je n 2 J E 

by definition (|1.36p and using the fact that W and Wk have the same minimum and mini- 
mizers (see the remark after (jl.34p ). We deduce 

F n {v n ) <n + a. 

Together with we find 1 ^ ^e" n 2 (r?+a) . Taking logarithms, we find 

logK^ > logn! + nlogJ 2 — nlogn — — n/3(rj + a). 

From Stirling's formula, logn! > nlogn — Cn and we deduce (|4.4p . with C v = — \og5 2 + C. 
Note that the dependence on rf comes from 5. □ 

5 Consequences for fluctuations: proof of Theorem [3] 

In this section, we prove Theorem [3] The first step is to find, via the method first introduced 
in [SS3] and tools from [SSI, ST] how F n and W control the discrepancy between u n and the 
measure n/xo- 

Proposition 5.1. Let v = Yli=i an d g v be as in Definition VJ.IA There exists a universal 
constant Rq > such that for any R > Rq, and any x' = \/tixq EM 2 , we have 

(5.1) f d 9u > cD(x' , Rf min ( 1, \E^±M) - CR 2 , 

Jb 2R (x' ) \ r j 

where c > and C depend only on V , and where D was defined in (fL38l) . 

4 In fact Ro could be any positive constant, and then c, C would depend on Ro as well, but this requires to 
adjust p accordingly and we omit for simplicity to prove this fact. 
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Proof. From Proposition 13.11 we have, defining j v as in Definition [3? 

(5-2) g v > -C + - A \j v \ 2 l 

where B p is a set of disjoint closed balls covering Supp(z/), and the sum of the radii of the 
balls in B p intersecting any given ball of radius 1 is bounded by p and we may take p < |. 

We distinguish two cases. Either D(x' ,R) > or D(x' , R) < 0. Let us start with the first 
case. We must introduce a modified distance function to x' , as follows. For any x, f(x) is the 
infimum over the set of curves 7 joining a point in Br(x' ) to x of the length of 7 \ B p . This 
is also the distance to x' Q for the degenerate metric which is Euclidean outside Br(x'q) U B p 
and vanishes on Br(x' ) U B„. We claim the following: 

Lemma 5.2. If \ X Xq I > R + 2 then 

|X ~ ^ - - < f{x) <\x- x' \ - R. 



Proof. The upper bound is obvious so we turn to the lower bound. Let j(t) be a continuous 
curve joining x = 7(0) to Br(x' ). Let us build by induction a sequence 



K 

JU *Aj ■ ■ • • ^ JLr 

with x k+1 defined as follows: let be the smallest t > such that j(t) ^ B(x k , 1) n 7 or 
j(t) E Br(x' ). This procedure terminates after a finite number of steps at x K G dBji(x' ). 
By triangle inequality we have 

(5.3) |x - x' | < \x k+1 - x k \ + - x' | <K + R. 

k=0 

On the other hand, by property of B p , for any < k < K — lwe have 

*(7[*fc,*fc+i] \6 P )>k fe+1 -x fc |-2p. 

Summing this over k and using (|5.3p . we find 

% \ P ) > if - 1 - 2ETp > |a? - x' | - i? - 1 - 2p (\x - x' \ - R) . 

Taking the infimum over all curves 7 we deduce that 

f(x) > {\x - x' | -R)(l- 2p) - 1. 

Since by assumption \x — x' \ — R > 2, we obtain f{x) > (\x — x' \ — R) (1 — 2p — i) and the 
result follows since p < 1/8. □ 

We proceed with the proof of Proposition 15. II Since / is Lipschitz with constant 1, almost 
every t is a regular value of /. For such a t the curve 7$ := {/ = t} is Lipschitz and does not 
intersect /3p, since V/ = there. Moreover, restating Lemma 15.21 we have 

(5.4) f{x) < t x € B fl+ «(so) U 5^+2(4), 

thus 74 C -B 2 J?(x' ) if R + it < 2R, i.e. if t < R/4, and R + 2 < 2.R, i.e. if i? > 2. It follows 
from (|5.2p and the coarea formula that if > 2 then 

(5.5) / dg v > -CR 2 + \ [ m ( [ \j u \ 2 ) dt. 

JB 2R (x' ) 4 Jt=0 \J-yt J 
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We proceed by estimating the innermost integral on the right-hand side. If t > 1/2 then 
Br+2 C Br+m and using (|5.4|) we find, using Definition 13.21 and writing D instead of D(x' , R) 
in the course of this proof, 



Ju ■ t 

it Hf<t} 



j curij'j, > I B_r_(x ) ) — 2irnfio ( Br+u{xq) ) 
J{f<t} \ ^ J \ V* / 

> 2-kD -C((R + At) 2 - R 2 ) 



where we have used (|1.15p . The right-hand side is bounded below by irD if R+4t < V R 2 + cD 
and c is small enough. Thus, if 



(5.6) l/2<T:=|^l + c^-lj, 

it follows using Cauchy-Schwarz's inequality that for every t £ (1/2, T) we have 

n 2 



(5.7) / \j,\ Z >7T 



Inserting into (|5.5p while reducing integration to 1/2 < t < min(T, -R/4) we obtain 



r r 
(5.8) / d 5iy > -CR 2 + ir 2 D 2 I 

JBon(x' n ) Jt- 



mm(T,R/4) j 

: — rdt. 



B 2R {x' ) Jt=l/2 nt\ 

The evaluation of the last integral is complicated by the fact that 7t is the level set for / 
rather than the usual circle, however the result will be comparable to the one we would get 
if we had |7t| = 2ir(R + t), this is proven as follows: From Lemma 15.21 we have for every 
t € [l/2,min(T,jR/4)] that j t C {x : < \x - x' \ - R < min(4T,i?)}. From the coarea 
formula and the fact that |V/| < 1 it follows that 

/ ht\dt< \{x:R< \x-x' \ < R + mm(4T,R)}\ <CRmm(4T,R). 

Then, using the convexity of x i— > 1/x and Jensen's inequality in (|5.8[) we obtain for some 
c > that 

,min(T,R/4) 1 ( m in(T, R/4) - i) 2 

'l l7t[ ^ ~ C i?min(4T, i?) ' 

Inserting into ()5.8[) we obtain assuming (|5.6p and 

(5.9) Kmin(T,i?/4) 
that 

(5.10) / dg v > -CR 2 + f mm(T, R/4) - )\ . 
Jb 2R ( x ' ) U V 2 / 

One may check that (|5.6p . ([5.9p are satisfied if R > 4 and -D > Co-R for a large enough 
Cq > 0. Then it is not difficult to deduce (|5.ip from (|5.10p by distinguishing the cases T < R/4 
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and T > R/A, i.e. D < C X R 2 and D > C X R 2 for a well chosen C x . Finally, if D < C R then 
(|5.ip is trivially satisfied, if C is chosen large enough. 

Let us turn to the case D(x' ,R) < 0, which implies \D(x' Q ,R)\ < rifio(B R /^(xo)) < 
irR 2 m. 

In this case we define f(x) to be the distance of x to the complement of Bjt(x' Q ) with 
respect to the metric which is Euclidean on Br(x' ) \B p . As in Lemma 15.21 and assuming 
\x — x' Q \ < R — 2 we have 

- A ° A < /(*) <R-\x- x' \. 

and as above, if R > 2 and for almost every 1/2 < t < R/A the curve Jt = {f = is a 
Lipschitz curve which does not intersect B p and {/ < t] C Br(x' ) \ BpL-itix'o)- It follows 
that, writing as before D for D(xq,R) 



(5.11) / j u -T= I curl> = / curlj,, - / curl> < 

Jit J{f<t} JBr{x' q ) J B(x' ,R)\{f>t} 

2-kD + 27rn/x (5 JL \ Br-u ) < 2ttD + C (R 2 - (R — it) 2 ) . 

The proof then proceeds as in the first case by using Cauchy-Schwarz's inequality and inte- 
grating with respect to t £ [1/2, min(T, R/A)], where 



which ensures that the right-hand side in (|5.1ip is bounded above by irD. Note that D is 
nonpositive, but bounded below by —CR 2 hence if c > is small enough the quantity inside 
the square root above is positive. □ 

We now proceed to the proof of Theorem [31 starting with (jl.42p . If R > Rq and 
\D(x'q,R)\ > rjR 2 then from Proposition 15.11 and using the fact — from Proposition 13.11 
- that g u is positive outside U'^ =1 B(x' i , C) and that g v > —C everywhere, we deduce from 
([377)1 and dSU) that 

(5.12) F n (v) > i (-CR 2 + cmm(r] 2 ,r] 3 )R 4 ) + 2 J (dv. 

Inserting into (jl.23p we find 

FP(\D(x' ,R)\ >t]R 2 ) <-^exp(Cf3R 2 -c/3mm(r] 2 ,r] 3 )R A ) j e' n ^ ^ dv dx x . . . dx n . 

Then, using the lower bound (|4.4p and Lemma 13.51 we deduce that if /3 > /?o > and n is 
large enough depending on (3q then 

logpP(\D(x' ,R)\ >7]R 2 ) <-c(3mm(ri 2 ,ri 3 )R 4 + CPR 2 + CnP + Cn, 

where c, C > depend only on V. Thus (|1.42|) is established. 
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We next prove (jl.44p . By Fubini's theorem, and using again the facts that g u is positive 
outside U^ =1 B(x' i ,C) and > — C everywhere we have 

/ dg v > [ l-f dg v ) dx' - C\U'\ - Cn. 

JR 2 JU' \Jb(x',2R) J 

Combining with Proposition 15.11 it follows that 

J dg v > -C{\U'\ + n ) + ^J -CR 2 + cD(x', R) 2 min (l, dx'. 
i.e., changing the constants if necessary, 

(5.13) dg v > -C(\U'\ + n) + J D(x', R) 2 min (l, dx. 

It follows, using as above (|1.23]) . (|3.T|) . f|4.4j) and Lemma [331 and since \U'\ = n\U\, that 

logP^ (J min (l, |£>( ^^h dx > 77) < -cnf3n + Cn/3 (\U\ + 1) + Cn, 

where c, C > depend only on V, where (3 > (3q > and where n > no(A))- 

We next turn to ()1.43p . Arguing as above, from (|3.7p we have -F n (^) > — C + 2 J* C^dv. 
Splitting 2 f Qdv as f Qdv + J C^S inserting into (|1.23|) and using (|4.4|) we are led to 



where C depends only on V. Then, using Lemma 1331 we deduce (|1.43p . 

There remains to prove (|1.45|) . Let v = Y^1=i an d X be a nonnegative function such 
that x = 1 011 U ■= E' U (u™ =1 .B(a^, ^)) U -Br and ||x||oo> || Vx||oo < 1 5 compactly supported 
on J7 = {x : d(x, U) < 1}. We have \U\ < C(n + R 2 ), where C depends only on E, i.e. on V. 
Then Corollary 1.2 in [ST] asserts that, for any q < 2, 



(5-14) Wy/XJuWl* <C q (n + R 2 y«-i (W{j u , X ) + n)* . 

where we use that by construction, #{p G A\B(p, i) n {0 < X < 1} 7^ = 0. This same fact 
implies that v' = in the support of 1 — hence 

W{j v ,l- X ) = \ [ (i - x)\j,\ 2 > 0. 



In particular W(j u , X ) < > x) + W0V> 1 — x) = 1r 2 )- It then follows from (|5.14p 

and the fact that F n (v) = ^W(j u , 1 R2 ) + 2j(dv> (cf. (TOll ). that 



(5.15) y IjVI 9 < C^n + i? 2 ) 1 -^! \F n (v)-2J (dv + 1 

But, reasoning as above after (|5,12p . the probability that F n {y) — 2 § Qdv > n is bounded 
above for any (3 > (3q > an n large enough depending on (3q by exp(— in/3?] + Cn((3 + 1)), 
where C depends on V only. In view of (|5,15p . 

K(j B \j v \ q >C q n( V + l)i + 'I <exp(-in/3r ? + Cn(/3 + l)). 
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Rescaling we have Jq \VH \ q = 1 fQi\ju\ 9 , where H = — 2ttA 1 (v — nno), while 
nno\\ w -i, q{n) < C\\VH\\ Lq ^y Therefore 



>2> 1 ' 



P " (( 1 + ^") 2 9 \\ u - n to\\w-i«(B R/ ^)>C q n*(l+ri)*J <exp(-in/3 ?? + Cn(/3 + l)). 
After a slight rewriting, this concludes the proof of Theorem [3j 



6 Lower bounds via the ergodic theorem and conclusions 

6.1 Abstract result via the ergodic theorem 

In this section, we present the ergodic framework introduced in [SS1| for obtaining "lower 
bounds for 2-scale energies". We cannot directly use the result there because it is written 
for a uniform "macroscopic environment", which would correspond to the case where mo(x) 
is constant on its support (as in the circular law). To account for the possibility of varying 
environment or weight at the macroscopic, we can however adapt Theorem 3 of [SS 1J and 
easily prove the following variant: 

Let X denote a Polish metric space, when we speak of measurable functions on X we will 
always mean Borel-measurable. We assume there is a d-parameter group of transformations 
9\ acting continuously on X. More precisely we require that 

- For all u G X and A, // € M. d , 6\{6^u) = 6\ +ll u, 9qu = u. 

- The map (A, u) >— > 9\u is continuous with respect to each variable (hence measurable 
with respect to both). 

Typically we think of X as a space of functions defined on M. d and as the action of trans- 
lations, i.e. 6\u{x) = u(x + A). Then we consider the following d-parameter group of trans- 
formations T| acting continuously on M, d x X by T^(x,u) = (x + e\,6\u). We also define 
T x (x,u) = (x,6\u). 

For a probability measure P on R d x X we say that P is translation-invariant if it is 
invariant under the action T, and we say it is T^^-invariant if for every function A(x) of class 
C , it is invariant under the mapping (x,u) i— > (x,9x( x ) u )- Note that T^^-invariant implies 
translation-invariant. 

Let G denote a bounded compact set in M. d such that 
(6.D |G|>0, lin J ( ° + g AG| = 0, 

for every x E M. 2 . We let {f £ } £ and / be measurable nonnegative functions on G x X, and 
assume that for any family {(x £ ,u £ )} e such that 

\/R > 0, limsup / f £ (T{(x £ , u £ )) d\ < +oo 
£->o Jb r 

the following holds. 

1. (Coercivity) {(x £ ,u £ )} £ admits a convergent subsequence (note that {x £ } £ subsequen- 
tially converges since G is compact). 
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2. (r-liminf) If {(x £ ,u £ )} £ converges to (x,u) then 

liminf f £ (x e ,u £ ) > f(x,u). 

e— >0 

Then we consider an increasing family of bounded open sets {Ui?}i{>o such that 

(6.2) (i) {XJr}r >0 is a Vitali family, (ii) lim K A + ^fl) AU fl | = Q 

i?-s>+oo |Ur| 

for any A G M. d , where Vitali means (see |Rij ) that the intersection of the closures is {0}, that 
R i y |Ur| is left continuous, and that |Ur — Ur| < C\XJr\. 
We have 

Theorem 6. Let G, X , {9\}\, and {f £ } £ , f be &s above. Let 

F £ {u) = i f e (x,9*u)dx. 
Jg 

Assume that {F e (u £ )} e is bounded. Let P £ be the probability on G x X which is the image of 
the normalized Lebesgue measure on G under the map x h-> (x,9s.u £ ). Then P £ converges to a 
Borel probability measure P onGxX whose first marginal is the normalized Lebesgue measure 
on G, which is Txm -invariant, such that P-a.e. (x,u) is of the form lim £ ^o(x £ ,9s^u £ ) and 
such that 

(6.3) liminf F £ (u £ ) > / f(x,u)dP(x,u). 

e^O J 

Moreover, 

(6.4) [ f(x,u)dP(x,u) = E p ( lim / f(x, x u) dX ) , 

J \R->+oaJ l j R J 

where E p denotes the expectation under the probability P. 

Proof. We only sketch the proof as it is very similar to [SSI j . The following points have to 
be checked: 

1. P £ is tight hence has a limit P. This follows from the coercivity property of f £ as in 

2. P is -invariant. Let $ be bounded and continuous, and let P\ be the push-forward 
of P by (x,u) h-» (x,9x( x )u). Then from the definition of P\, P, P £ , we have 

<S>(x,u)dP x (x,u) = j §(x,9 xfyX) u)dP(x,u) = lim J $>{x,9 x{x) u) dP £ (x,u) = 

f f *((! + e\)- 1 (y),9vu £ ) 

lim + 9i._i_\( T \u E ) dx = lim + — — — — —. . E . . ... dy, 

e-*oJ G l ' c +x{x) £> £ ^oJ {I+eX)iG) \det(L + eD\((L + e\)-i(y))\ 

where the last equality follows by the change of variables y = (J + e\)(x). Using the 
boundedness of <&, the C 1 character of A, the compactness of G and (|6.ip . we may 
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replace (I + e\)(G) by G and the denominator by 1 in the last integral and we find, 
using the definition of P £ 



(6.5) J <5>(x,u)dP x (x,u) = lim J $((/ + eA) _1 (x), u) dP £ (x, it). 

Since {P E }s is tight, for any S > there exists Kg such that P £ (Kg c ) < 5 for every e. 
Then by uniform continuity of $ on Kg the map (x,u) i— > < &((7 + eA) _1 (x),u) converges 
uniformly on to (sc, u) h- >• and thus 

lim / <J>((I + eA) -1 (x),u)dP e (£,ii) = lim / $(x,u) dP e (x,u). 

Since this is true for any 5 > 0, and using the boundedness of $ we get 

lim / ^((1 + eXy 1 (x),u) dP £ (x,u) = lim / $(x, it) oLP £ (x, u) = / $(x,u) dP(x,u), 
e-X)J e-*oJ J 

by definition of P. Thus in view of (|6.5p we have P\ = P and P is thus Tx^-invariant. 

3. liminf^o F £ (u £ ) > J fdP. This follows from Lemma 2.2 of |SS1| . since F £ (u £ ) = 
// '//'• 

To conclude, as in |SSlj . the fact that P is T^^-invariant (which implies T\-invariant) and 
Wiener's multiparametric ergodic theorem (see e.g. |Bej ) imply that 

f f(x,u)dP(x,u) = E p f lim / f(T x (x,u))d\ ) = E p ( lim / /(x, # A u) dA ) . 

□ 



6.2 Proof of Theorem [2], part A 

The proof follows essentially [SSI] . Proposition 4.1 and below. Let {f n }n and P Vn be as in the 
statement of Theorem [21 We need to prove that any subsequence of {P Un }n has a convergent 
subsequence and that the limit P is a X^^-invariant probability measure such that P-almost 
every (x, j) is such that j S A mo ^ and (|1.39j) holds. Note that the fact that the first marginal 
of P is dx\E/\E\ follows from the fact that, by definition, this is true of P Un . 

We thus take a subsequence of {P Un } (which we don't relabel). We may assume that it 
has a subsequence, denoted u n , which satisfies F n {D n ) < C, otherwise there is nothing to 
prove. This implies that D n is of the form Y^h=i ^ « ■ ^ e ^ ^ n denote the current and g n the 
measures associated to i> n as in Definition 13.21 As usual, v' n = Y!2=i ^yjnxi „ ■ Another useful 
consequence of F n {v n ) < C is that, using ()2.ip . we have w n (xi^ n ) = n 2 I(fio) + 0(n log n), 
which in turn implies (see |ST| ) that 

(6.6) -v n -)• /j, . 

n 
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Step 1: We set up the framework of Section 16.11 

We will use integers n instead of e to label sequences, and the correspondence will be e = 
l/x/n. We let G = E and X = M + x Lf oc (M 2 , M 2 ) x M, where p G (1, 2), where M+ denotes 
the set of positive Radon measures on M? and A4 the set of those which are bounded below 
by the constant — Cy := — C(||mo||oo + 1) of Proposition ^. 1| both equipped with the topology 
of weak convergence. 

For A £K 2 and abusing notation we let 8 X denote both the translation x — > x + A and the 
action 

Accordingly the action T n is defined for A E M 2 by 

TZ(x,v,j,g) = U + -j=,9 x #v,jo6 x ,6 x #g) . 



Then we let x be a smooth cut-off function with integral 1 and support in -6(0, 1) and define 

fa „, , , \ \ f x(y)dg(y) X(vj,g)=0> iia .(p f n ,j n ,gn), 

(6.7) f n {x,v,j,g) = <j vr J R 2 

[ +oo otherwise. 
Finally we let, in agreement with Section f6.1| 

(6.8) F n (u, j, g) = f n (x, 9 x ^(u, j, gfj dx. 
We have the following relation between F n and F n , as n — > +oo: 



(6.9) F n (u, j,g) is 



^ ]h\ { F n{v n ) - 2f(dv n ) + o(l) if (v,j,g) = (v' n ,Jn,g n ) 
= +oo otherwise. 



Indeed it is obvious from (|6.7p that if (v,j,g) (v' n ,3n,gn) then ~F n (v,j,g) = +oo. On the 
other hand, if (f,j,g) = (v' n ,jn,g n ), then from the definition of the image measure 6 x #g n , 

~F n (v,j,g) = ^J-^J x(y ~ Xy/n) dg n {y) dx = J x*^E'dg n . 

Since x*\e> is bounded above by 1 and is equal to 1 on U := {x' : dist(a/, M 2 \ E') > 1} we 
deduce that 

(6 ,o) ,p. ( ,, A ,) < tffj-g^ = "-w+gvn 



F n {D n )-2 J Qdu^j + 



It— t I \ — lb lb / ~ I > ""ft I i It— ■ i 

\E\ \ J J n\E\ 

Then we note that from (|3.2p in Proposition 13.11 the measure g~ is supported in the union of 
balls B(x', C) for x' G Supp(z^), and bounded above by a constant. Thus g~(U c ) is bounded 
by a constant times the number of balls intersecting U c , hence by Cv' n \x' : dist(x', U c ) < C}. 
From (|6.6p this is equal to Cnno{x : dist(x, dE) < C/y/n} + o(n), which from the assumption 
(|1.13p on E is o(n). Plugging this into ()6.10p proves ()6.9p . 
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Step 2: We check the hypotheses in Section 16.11 

We must now check the Gamma-liminf and. coercivity properties of {fn}n* 
Lemma 6.1. Assume that {(x n , v n , j n , g n )} n converges to (x,v,j,g). Then 

liminf f n (x n ,v n , j n ,9n) > f{x,v,j,g) := - [ xdg. 

n ir J 

Proof. We may assume that the left-hand side is finite, in which case fn(x n ,v n ,j n ,g n ) = 
J xdg n for every large enough n, from which the result follows by passing to the limit. □ 

Lemma 6.2. Assume that for any R > we have 

(6.11) limsup [ i n ( x n + -^=,6 x(vn,jn,gn)) d\ < +oo. 

n->+oo J B R V V n / 

Then a subsequence of {(x n , v n , j n ,gn)}n converges to some (x,u,j,g) G E x X. 

Proof. Assume (|6. 1 1[) . Then the integrand there is bounded for a.e. A and from the definition 
(16. 7p we deduce that 

0\(Vn,jn,9n) = 9y/jix n +\(v'nJn,9n) 
and then that (u n ,j n ,g n ) = 9 v ^ Xn (v' n ,j n ,gn)- Thus (|6.1ip gives, in view of (|6.7p . that for 
every R > there exists Cr > such that for any n 

j B j X(y~ Vnx n - A) dg n (y) dX = J x* 1 B R (^Ex n ) dg n < C R . 

This and the fact that g n is bounded below implies that g n (Bji(y/nx n )) is bounded in- 
dependently of n and then, using (|3.4|) . that the same is true of P' n (Bn(s/nx n )). In other 
words {u n = 9^ Xn v' n } n is a locally bounded sequence of (positive) measures hence converges 
weakly after taking a subsequence, and the same is true of {g n = d^/n Xn g n }n- On the other 
hand {x n } n is a sequence in the compact set E hence converges modulo a subsequence. 

It remains to study the convergence of {j n = jn°9^/E Xn +\}n- From (|3.3|) in Proposition ^. II 
and the local boundedness of {v n }n we get that W(j n ,x* ^B R (^Ex n )) = W(j n ,x * 1b r ) is 
bounded independently of n for any R > and then, using Corollary 1.2 in [ST], that {j n }n 
is locally bounded in L^ C (R 2 , M 2 ), for any 1 < p < 2 hence a subsequence locally weakly 
converges in L^ c (IR 2 , M 2 ). Strong local convergence follows as in |SSlj using elliptic regularity 
with curlj n = v' n — \i§ and div j n = 0. □ 

Step 3: Conclusion 

From the previous steps, we may apply Theorem [6] in this setting and we deduce in view of 
(|6.9p that, letting Q n denote the push- forward of the normalized Lebesgue measure on E by 
the map x i-> (x, 9^ ix {v' n ,j n ,g n )) 1 and Q = lim n Q n , 

(6.12) liminf ^F n (v n ) - 2 J (dv n ^J > liminf F n (i>' n ,j n ,g n ) > 

J {^j Xdg^j dQ(x,u,j,g) = J^hm^^ J x(y ~ A) dg(y) d\dQ(x,u,j, g) = 

Now we use the fact that for Q-almost every (x,is,j,g): 
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i) There exists a sequence {x n } n in E such that (x,u,j,g) = lim n (x n , 0^ lXn (u' n , j n , g n )). 

ii) As a consequence of the above — f x*^-B R dg converges to a finite limit as R — > +00. 

The first point implies, since curlj n = v' n — and div j n = 0, that curlj = v — mo(x) and 
div j = 0. The second point implies in particular using (|3.4p that v(Br) < CR 2 , proving 
that (j» G AnoW- 

Moreover the second point implies that for any C > we have g{B R+ c \ Br-c) = °{B 2 ) 
as R — > +00, and thus from point i) above 

lim lim ^g n {B R+c {\/nx n )) \ B r ^ c {Vnx n )) = 0. 

Using ()3.4p we deduce that 

lim lim -^V>' n (B R+c {Vnx n )) \ B r _ c (Vnx n )) = 



and then from (|3.3[) . 



lim lim — - 

R-^+oc n-5-+oo ij 2 



Thus, using Lemma 4.8 in |SS1] in the limit with respect to n we deduce 



0. 



0. 



lim W ti> X* 1 B R )~ [ x*lB R dg 
Together with (|6.12p this yields, by definition of W, 

(6.13) liminf ±- f F n (i/ n ) - 2 j Qdu^j > j W(j)dQ(x,u,j,g) 

and, we recall, Q-a.e. (J, u) £ A mo r x y 

Now we let P n (resp. P) be the marginal of Q n (resp. Q) with respect to the variables 
(x,j). Then the first marginal of P is the normalized Lebesgue measure on E and P-a.e. we 
have j £ A mo ( x ), in particular 

W(j) > min W = mo(x) ( minVF — — logmo(x) j . 

Integrating with respect to P and noting that since only x appears on the right-hand side 
we may replace P by its first marginal there we find, in view of (|1.36p that the lower bound 
(fj~39l) holds. 



6.3 Proof of Theorem [2], completed 

As mentioned above, Part B of the theorem is a direct consequence of Proposition 14.11 see 
Corollary IP 

Part C follows from the comparison of Parts A and B and the fact that W and Wk 
have the same minimizers and minimal value. This implies that, for minimizers, the chains 
of inequalities (|1.39p and (jl.40p are in fact equalities and that J W dP must be minimized 
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hence equal to a. Since (jl.39p follows from (|6.13p . the inequality (|6.13p must also be an 
equality and moreover we must have lim Tl _ ) . 00 j (,dv n = 0, which in view of (|1.17|) . implies 
that lim]>]idist(xi,£) 2 = 0. 

From the fact that F n (v n ) = O(l), we deduce from Proposition EID and (|3.7|) . (arguing 
exactly as in the proof of Theorem [3]) that there exists C > such that for every x, R > 1, 
we have 

D(x, R) 2 min (l, < Cn. 

This completes the proof of Theorem [2j 



6.4 Deviations: proof of Theorems [5] and Theorem Q] 

We start with the upper bound on logP^. Let A n be a subset of <C n . We identify points in 
C n with measures v n of the form Y2i=i *W 
From (|1.23|) . we have 

F^(A n ) = \\ e-3^E*=l dx X ... dx n 
Kin J A n 



hence 



(6.14) 



\og¥ p n {A n ) \ogRP , 1 



n n 



+ - log f e -5^ F »(E? =1 dxi . . . dXn . 
■J An 



We deduce, writing F n (v) = F n (u) — 2 J C^dv^ that 

(6.15) lQgP ' (A) < + I log (e-¥^ An ? n f e -P»fC*n dxi _ dx \ . 

n n n \ J An ) 

Let v n such that F n (y n ) < inf^ n P n + 1/n. Then from (|1.39p in Theorem [2] we have, using 
the notations there, liminfn-^oo F n (y n ) > ^ J W(j) dP(x, j) where P = \\m n P Un . Since 
P Vn G i(^4 n ) by definition we have P G since = n n> o { J m>n i(A m ) . We may thus write 

\E\ f 

(6.16) liminf F n {v n ) > J— I- inf W(j)dP(x,j). 
Inserting into (|6.15p we are led to 

logP£(A n ) ^ /3|P| . . f nr ,.. JT1 , logi^ 



+ -log( [ e~ f)n ^ dVn dxi...dx^\ +o(l) 



thus in view of Lemma 13.51 and (|4.4p , we have established (|1.47p . An immediate corollary of 
(|6.17|) . choosing A n to be the full space and using inf — f W(j) dP{j) = a and Lemma 13.51 is 
that 

log Kn 0a , , „, 
(6.18) limsup & <-— + log|P|. 
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We next turn to the lower bound. Fix rj > 0. Given A, let P 6 A be such that 



(6.19) 



f W(j)dP(x,j) < inf f W(j)dP(j) + %. 
J PeAJ 2 



Since P £ A, if r\ is chosen small enough (which we assume) then B(P,2r]) C A, where the 
ball is for a distance metrizing weak convergence as in Proposition 14. 11 

We then apply Proposition 14.11 to P and r]. We find 5 > and for any n large enough a 
set A n such that \A n \ > nl^S 2 /n) n and, rewriting (|4.1|) with (|2.1|) . 



(6.20) 



\E\ f 

limsupsupF n < — / W K {j)dP(j) +rj. 
n->oc A n T i 



Moreover, for every (yi,...,y n ) G ^4„ and letting {i/„, = 5^ H h (5 yn }„, there exists {j„}„, 

in Lf oc (M. 2 , M 2 ) such that curlj n = 27r(u' n — mo') and such that the image P n of cLx\e/\E\ by 
the map x \-t (x, j n (\/nx + •)) is such that 

(6.21) limsupdist(P n ,P) < rj. 

In particular (|1.49p holds. Moreover, inserting ()6.20p and (|6. 19|) into (jl.23fl . we find that 

A r . 



logP^(^ n ) > JogK^ _ p\E[ 
n ~ n 2n 



inf / W{j)dP{j)-\pr ] + -log 



+ 0(1). 



pga- 

On the other hand, using \A n \ > n\(ir5 2 /n) n and Stirling's formula, we have log|A n | > 
2nlog5 - Cn. Combining with ([6TTH]) . (TQ5D follows, with C v = -2 log 5 + C + log 
Theorem [1] immediately follows by combining (|6.18[) . (|4.4p and ()1.24D . 



7 Proof of Proposition 14.11 

The construction consists of the following. We are given e > 0, which is the error we can 
afford. First we select a finite set of currents J\, . . . , Jn (N will depend on e) which will 
represent the probability P(x,j) with respect to its j dependence, and whose renormalized 
energies are well-controlled. Since P is T\^ x ) -invariant, we need it to be well-approximated by 
measures supported on the orbits of the Jj's under translations. Secondly, we work in blown- 
up coordinates and split the region E' (whose diameter is order y/n) into many rectangles K 
with centers xk and sidelengths of order R large enough. Even though we choose R to be 
large, it will still be very small compared to the size of E' , as n — > oo, so that the Diracs 
at xx/y/n approximate P(x,j) with respect to its x dependence. On each rectangle K, the 
weight mo' is temporarily replaced by its average mjf. Then we split each rectangle K into 
q 2 identical rectangles, with sidelengths of order 2R = R/q, where both R and q will be 
sufficiently large. We then select the proportion of the rectangles that corresponds to the 
weight that the orbit of each Jj carries in the approximation of P. In these rectangles we 
paste a (translated) copy of Jj at the scale tuk and suitably modified near the boundary 
according to a construction of [SSI] (Proposition 17.41 below) so that its tangential component 
on the boundary is (this can be done while inducing only an error e on W). In the few 
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rectangles that may remain unfilled, we paste a copy of an arbitrary Jo whose renormalized 
energy is finite. We perform the construction above provided we are far enough from dE' . The 
layer near the boundary must be treated separately, and there again an arbitrary (translated 
and rescaled) current can be pasted. Finally, we add a vector field to correct the discrepancy 
between mx and uiq in each of the rectangles. 

To conclude the proof of Proposition l4.ll we collect all of the estimates on the constructed 
vector field to show that its energy w n is bounded above in terms of J W dP and that the 
probability measures associated to the construction have remained close to P. 



7.1 Estimates on distances between probabilities 

First we choose distances which metrize the topologies of Lf oc (M 2 ,IR 2 ) and B(X), the set of 
finite Borel measures onI = £x Lf oc (M 2 , M 2 ). For j x , j 2 G Lf oc (M 2 ,M 2 ) we let 

II jl ~ h\\LP(B{0,k)) 



dp(jx,32) = J2 2 



k=l 



1 + II Jl - J2\\ LP (B(0,k)) 



and on X we use the product of the Euclidean distance on E and d p , which we denote dx- 
On B{X) we define a distance by choosing a sequence of bounded continuous functions {(fk}k 
which is dense in Cb(X) and we let, for any fj,i, [A2 6 

At,, ,^_X^o-fc l(yfc^i-M2)| 

where we have used the notation ((p, n) = J ip dfi. 
We have the following general facts. 

Lemma 7.1. For any e > there exists r/Q > such that if P,Q £ S(X) and \\P — Q\\ < 7/0, 
then d(P,Q) < e. Here \\P — Q\\ denotes the total variation of the signed measure P — Q, i.e. 
the supremum of (<p,P — Q) over measurable functions ip such that \<p>\ < 1. 

In particular, if P = Y^Li a i$xi and Q = Z^iA$c< with £V \cti - /%| < r) , then 
d B (P,Q) <e. 

Lemma 7.2. Let K C X be compact. For any e > there exists n\ > such that if 
x 6 K,y £ X and dx(x,y) < rji then ds(5 x ,5 y ) < e. 

Lemma 7.3. Let < e < 1. If fx is a probability measure on a set A and f,g:A—*X are 
measurable and such that dj3(Sf^,8 g ^) < e for every x £ A, then 

d B (f*n,g#v) <Ce(|log e| + l). 

Proof. Take any bounded continuous function ip^ defining the distance on B{X). Then if 
&g(x)) < e f° r an Y x £ X we have in particular 

Wk{f{x)) -ipk{g{x))\ <2 k £ 



1 + Wk{f{x)) - (p k (g(x))\ 

It follows that 

00 

d s (/# M ,ff # M) <^2- fc min(e2 fc ,l) < e ([log 2 e ] + 1) + £ 2' k < Ce(|log e\ + 1) 

k fc=[log 2 e]+l 



□ 
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7.2 Preliminary results 



In what follows E' = \/nE, mo'(x) = mo(x/y / n): we work in blown-up coordinates. We 
consider a probability measure P on Ex Lf oc (M 2 , M 2 ) which is as in the proposition. We let P 
be the probability measure on E x Ai which is the image of P under (x,j) \— > (x, o'i/ mQ ( x \j)i 
so that 



(7-1) P = J S x 5 ai/mo{x)j dP(x,j), P = J 6 X ® 5 amo{x)j dP(x,j) 

It is easy to check that since P is T^^-invariant, P is as well, and in particular it is translation- 
invariant. 

The construction is based on the following statement which is a rewriting of Proposition 4.2 
in |SSlj and the remark following it: 

Proposition 7.4. Let Kr = [— R, R] 2 , let {xr}r satisfy (|1.29|) and let Wk be the renormal- 
ized energy relative to the family {Kr}r. 

Let G C A\ be such that there exists C > such that for any j G G we have 

(7.2, ™>i.niH 

for the associated v 's, and 

(7.3) lim ^fl=W K (j), 

R^+co \Kr\ 

where the convergence is uniform w.r.t. j € G. 

Then for any e > there exists Rq > 0, 772 > such that if R > Rq and L is a rectangle 
centered at whose sidelengths belong to [2R,2R(1 + 7/2)] and such that \L\ £ N, then for 
every j E G there exists a jl E L^ oc (IR 2 , M?) such that the following hold 

i) 3L = in L c , 

ii) There is a discrete subset A C L such that 

curl j L = 2tt \y~^ S p - 1 L 

\peA 

Hi) Ifd(x,L c ) > Ri then j'l(x) = j{x) 



iv) 



(7.4) W(jL L > 1l) <W K (j) + e. 



The next lemma explains how to partition E into rectangles. 
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Lemma 7.5. There exists a constant Co > such that, given any R > 1 and q G N*, i/iere 
exists for any n G N* a collection fC n of closed rectangles in E' with disjoint interiors, whose 
sidelengths are between R = 2qR and R + CqR/R 2 , and which are such that 

{x <E E' : d(x, dE') <R}CE'\ \J K C {x G E' : d(x, dE') < C R} 

KeKn 

and, for all K G JC n , 

(7.5) / m ' G g 2 N. 

Proof. For each j 6 Z we let 

m j{t) = / / _ mo'(x, y) dy dx. 

J x=—oo J y=jR 

Then each strip {jR < y < (j + l)R} is cut into rectangles [Ujit(i+i)j] x \jR, (j + where 
toj = — oo and 

fi+lj = min{t > tij + : mj(tij) G g 2 N}. 

Since by assumption (|1.12j) we have mo'(x) G [m, fn] for any x G -E 1 ', it is not difficult to check 
that if such a rectangle is included in E' then 

- q 2 

Uj + R < ij+i j < tij + -R H =, 

mR 

and thus its sidelengths are between R and R+CR/R 2 since R/R 2 = Aq 2 /R. We let /C n be the 
set of rectangles of the form [%, t(i+i)j] x [j-R, + which are included in {x : d(x,dE') > 
R}. From the above, it follows that these rectangles in fact cover the set {x : d(x, dE') > CR} 
for some C > independent of R > 1, q. By construction each K G /C n is such that 

to ' = mj(i(j + iw) - mj(tij) G <? 2 N. 

□ 

The next lemma explains how to select a good subset of L^ C (IR 2 , R 2 ). 

Lemma 7.6. Let P be a translation invariant measure on X such that P-a.e. j G A\. Then 
for any e > 0, there exist subsets H £ C G £ in L^ oc (IR 2 , M?) which are compact and such that, 
for any R £ > 0, 

i) r]Q being given by Lemma \7.1\ we have 

(7.6) P(E x G £ c ) < min(r ?0 2 ,??oe), P(E x H c e ) < min(7 ?0 ,e). 

ii) For every j G H £ , there exists T(j) C K m R e such that 

(7.7) \T(j)\ < CR £ 2 m and A $ T(j) => 9 x j G G £ . 
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Hi) The convergence in the definition ofWxij) is uniform w.r.t. (x,j) G G £ and, writing 
curl j = 2tt(v — 1), 

(7.8) WkU) and are bounded uniformly w.r.t. (x,j) G G £ and R > 1. 

iv) We have 

(7.9) d B (P, P") < Ce(|log e| + 1), w/iere 

1 

/exff £ rn (x)\K Re \ J^^)K Re \r(j) 



P" = I - jzxnp i / . s x^S amQix)9tlj dndP{x,j). 



Moreover, there exists a partition of H e into U^P/* satisfying diam (H*) < 773, where r/3 is 

(7.10) j G PT £ , d p (j,f) <rj 3 , m G [m,m], /i G V^%\r(j) =^ d B (b~o m 6^,o~<j m 6^') < e; 
and i/iere exists for all i, Ji G P^ £ suc/i i/ioi 

(7.11) WW) < infWA: + e. 

Proof. 

S'iep P Choice of G £ . Since L^ c (M 2 ,]R 2 ) is Polish we can always find a compact set G, 
satisfying (|7.6p and P(G £ C ) < 770. Then from Lemma 17. 11 P\-G £ (the restriction of P to G, 
satisfies dg(P, PlG £ ) < e. 

From the translation invariance of P and for any A, we have P(E x 9\G £ ) > 1 — r/o 
and therefore d B (P, Pl.8\G £ ) < e. In view of (|7.ip . it follows that for any A £ I 2 we have 
||P — P\|| < r/o and then dg(P, P\) < e, where 



e 
e) 



Pa = / d~x ® ^ j dP(x,i) = / o" x ® £<? A<T , ^ dP(x,j). 

JExe x Ge JExGc 

Then using Lemma 17.31 we deduce that if j4 C M 2 is any measurable set of positive measure, 
then 

(7.12) d B (P,P')<Ce(\logs\ + l), where P' = I -f 5 X ® Sg^,^ d\dP(x, j). 

JexgJa 

Moreover, since P is T^^-invariant, choosing x to be a smooth positive function with 
integral 1 supported in B(0, 1), the ergodic theorem (as in Section [6TT1 or see again |Be| ) 
ensures that for P-almost every (x,j) the limit 



lim j±- [ W(j(X + -), X (X + -))dX 
-^+00 \K R \ J Kn 



r-^+oo \K R \ j k 

exists. Then 1r r * X is a family of functions which satisfies ()1.29[) with respect to the family 
of squares {Kr}r, and from the definition of the renormalized energy relative to {Kr}r we 
may rewrite the limit above as 

(7-13) W K (j)= lim -L-W(j,lK R *x)- 

R-^+oo \Kr\ 
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By Egoroff 's theorem we may choose the compact set G £ above to be such that, in addition 
to ([7.12p . the convergence in (|7.13|) is uniform on G £ . In fact, since Wk(J) < +00 and 
limsup R v{Kr)/R 2 < +00 for P-a.e. (x, j), where curlj = 2tt(v — 1), we may choose G £ such 
that dZS} holds. 

The arguments above show that the properties (|7.12|) . (|7.8[) can be satisfied for a compact 
set G e of measure arbitrarily close to 1. We choose G £ such that (|T.6j) holds. 

The next difficulty we have to face is that 6\j need not belong to G e if j does. 

Step 2: Choice of H £ . For j G G £ , let be the set of A's in \/WKr e such that 9\j g" G £ . 
Since, from (|7.6j) and the translation-invariance of P, for any A G I 2 we have P(Ex9\(G £ ) c ) < 
rjo 2 , it follows from Fubini's theorem that 

I \T(j)\dP(x,j) = [ P(Ex (9 x G £ ) c )d\ < 4mi? £ 2 min(r ?0 2 ,r ?0 e). 
Therefore, letting 

(7.14) H £ = {j G G £ : \T(j)\ < AmR^o}, 

we have that (|7.6p holds. 

Combining fj7.6f) and (|7.14p with Lemma 17.11 we deduce from f)7.12[) that ()7.9p holds, 
where we have used the fact that 9\o~ m j = o- m 9^^ x j to change the integration variable to 
H = v / m (a;)A in (|7TT2l . 

Step 3: Choice of J\, . . . , Jjy e ■ We use the fact that G £ is relatively compact in Lf oc (E?, M 2 ), 
Lemma [721 an d the fact that (m,j) 1— > o~ m j is continuous to find that there exists 7/4 > 
such that for any m G [m, fn] and any j G G £ it holds that 

(7-15) d p {j,j')<r} A =» (5 CTm j,^ m j') < e - 

Moreover, from the continuity of 1— > 9^j, there exists 773 > such that 

(7.16) j G G £ ,d p (j,f) < 773, \i G VrnK R£ =^ ^ (9^,9^') < r/ 4 . 

If j G -£f e and /j G then G G e hence applying (|7.15p . we get (|7.1U|) . 

Now we cover the relatively compact .£f e by a finite number of balls B\, ... , Bjv e of radius 
7/3/2, derive from it a partition of H £ by sets with diameter less than 7/3, by letting H\ = 
B-x n if e and 

HI+ 1 = B i+1 n H £ \ [B x U • • • U Bi) . 

we then have 

N e 

(7.17) H £ =\jHi, dmm (Hi) < m , 

where the union is disjoint. Then we may choose Jj G HI such that (|7.1ip holds. □ 
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7.3 Completing the construction 



Step 1: Choice of R £ . We apply Proposition 17.41 with G = G £ and y/mR, where m G [m, m]. 
The proposition yields rj2 > 0, R e > 1 such that for any j G G e and any m G [To) to] and any 
rectangle L centered at with sidelengths in [2y/mR £ , 2yfmR £ (\ + r/2)]) (|7.4p is satisfied for 
some jLt with it! replaced by y/mR £ . The reason for including yfm is that we will need to 
scale the construction to account for the varying weight mo(x). 

Since our rectangles will be obtained from Lemma 17.51 and we wish to use the approxi- 
mation by ji in them, we choose R £ large enough so that if m G [to, to] and L is a rectangle 
centered at zero with sidelengths in [2\/mR £ , 2y/mR £ (l + Cq/R £ 2 )] then 

(7.18) J^<V2, < Vo, i^^R^i-^) C {x : d(x, L c ) > \fmR e i } C K^^^, 

where Co is the constant in Lemma 17.51 C\ > 1 is to be determined later, and r/o is the 

constant in Lemma 17. 11 

1 ' 3 

If A G K^iRji-^) and since j = ji if d(x,L c ) > y/mR £ ' 1 , we deduce from (|7.18p that 

3 

#aJl = #AJ in -B(0, -y/mfie 4 ), so that from the definition of d p , taking R £ larger if necessary, 

(7.19) Vj G G £ ,m G [to,, to], A G K^ Re ^_ vo) , d p (8 x o- m j,6\o- m j L ) < —, 

where 771 comes from Lemma 17.21 applied on {o~ m j ■ TO G [to, to] , 7 G G £ }, i.e. is such that 

(7.20) m G [m, to], j G G £ ,f G Lf oc (R 2 , M 2 ) and d p (j,f) < Vl d B {6 0mj ,8 anj .) < e. 



Step 2: Choice of q £ and the rectangles. We choose an integer g e large enough so that 

(7.21) — ^ < r/o, ^ x max W^-((7 m Ji) < e 

m<m<m 

where Ci > 1 is to be determined later. We apply Lemma 17.51 with R £ , q £ and N £ to obtain 
for any n a collection tC n of rectangles (we omit to mention the e dependence) which cover 
most of E' , and we also apply Lemma \7~M We rewrite P" given by (|7.9p as 

(7.22) P"=Y1 L ~ ,1 K I \ , S X ^S e ,jd^dP(x,j). 



Now we claim that if n is large enough and x G Kj y/n, j G HI, /j, G y/mo(x)Kn £ \T(j), then 



(7-23) cfe ® ^ mo(x) 8^j,S XK ® Sa^e^ 



< 2s, 



where a;^ is the center of K/y/n and m^- is the average of too over K/y/m. Indeed, since 
too is C 1 we have |x — xk\ < C / y/n, \mo(x) — tuk] < C/y/n thus if n is large enough, since 
6nj G G £ we find 

d B (S x ® <W o0c) <U,<W ® <W K <? M j) < e. 
Moreover, since d p (j, Jj) < 773, we deduce from (|7.10p that 



dB ( S XK &> S^orftSnx ® 5 CTmjf e M j, ) < £, 
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which together with the previous estimate proves (|7.23p . 

Using (j7T2"3"j) together with Lemmas O EH and {777} , we deduce from (^|) and ([7722]) 
that d B {P,P"') < Cefllog e| + 1), where 



(7.24) 



where 



Ki<N s 



C( - r J \jm K K R 



Ki<N t 



(7.25) pi )jr = P[^=xHl 



v 



Step 3: Choice of subrectangles and current j n . We now replace pi t K in the definition (|7.24p 
by 



(7.26) 



|A| 



Qe 2 \E' 



jtU^k, where n i;K 



ql\E' 1 
|A| 



-Pi,K 



We have, since P(4- x Lf oc (M 2 , M 2 )) = |A|/|P'| 

A' 



(7.27) 
and 



x g £ 2 |P| , 



fe=l 



A 



I A 



IP' 



7]-»^,X - Pi,K 



< c 



xLf oc (M 2 ,R 2 ; 



|A| n»,K 



, 2 |P'| P 2 |P' 



Summing with respect to i and A, using the facts that Y^KeK l-^l < l-^'l' (E27J, and the 



fact that the cardinal of fC„ is 



\E'\ 



we find 



l<i<N s ,KeK„ 



I A 



P' 



He | 



We may always choose C\ large enough in (|7,18p and (|7.2ip so that the right-hand side is 
< rjQ. Then Lemma 1 7 . 1 1 implies that ds(P, P^) < Ce(|log s\ + 1) is still true after replacing 
Pi,K by ^jj^nijc in (|7.24j) . i.e. where 



(7.28) 



■ OuJi d^- 



Ki<N, 



Next, we divide each A £ /C n into a collection £^ of q £ 2 identical subrectangles in the 
obvious way and we partition Ck into collections Cxi, < i < N £ such that if k > 1 then 
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Ck.,% contains n^x subrectangles. This is clearly possible from (|7.27p . If the inequality is 
strict we put the extra subrectangles in Ck,o, there will be uq^k of them and then 



N e 

2 

£ 

k=0 



(7.29) y~] n k ,K = q, 

We rewrite ([T728]) as 

(7-30) P(4) = W\ £ t'xVt-W*- 

Ki<N e 



LeC 



K.i 



Now, for L G C>K,ii let L = yJrriK{L — x£), where denotes the center of L. From 
Lemma IT3| a rectangle K G JC n has sidelengths between 2q £ R £ and 2q £ R £ (l+Co/R £ 2 ). There- 
fore L is a rectangle centered at zero with sidelengths between 2^/m.KRe an d 2^/rnKR £ {\ + 
C /R £ 2 ), and (TTTHD holds. 

This, and the results of Lemma 17,6} allow us to apply Proposition 17.41 on L to any Jj, 
1 <i < N £ . Note that \L\ E N follows from the fact that 

\L\ = m K \L\ = -f m '^l = -Ay / mo' 
'if & <?e JK 

and (|7.5p . In this way, we define currents Ji t i which satisfy (|7.4p and ()7.19p . We claim that, 
as a consequence of the latter, we have 
(7.31) 



f = Jll on L => di 



(7 

\J^K~K Re rn K \K R£ \ J L J 

This goes as follows: (i) Using Lemma 17.11 and (|7.7p . (|7.18p . we find that integrating on 
y/rnKKn V0 \jl e \ r(Jj) instead of y'mjfK^ and L induces an error of Ce. (ii) From (|7.19|) . 
and (|7.20p applied to 9^3% and 6^f we have de(#0 j i5 ^j') < £ and thus in view of Lemma [T73l 
we may replace O^Ji by (Lj' in the integral with an error of Ce|log e| at most, (iii) Using 
(|7.18p . (|7.7p and Lemma 1 7 . 1 1 again . we may integrate back on ^rnxKn^ and L rather than 
on K(i_ V0 ^ Re \ r(Jj), with an additional error of Ce. this proves ()7.3ip . 

Combining flEEED with (f?30l) and d B (P, P< 4 )) < Cefllog e| + 1), using Lemma 17.31 we find 
d B (P,P {5) ) < Ce(|log e| + 1), where 
(7.32) 

\E'\ K ^ m K J L XK °™kW* * \E'\ ^ J L/VW ^ XK ^m K J i: L 

l<i<N e l<i<N e 



where the last equality follows by changing variables to A = /j,/y/rn K , and where denotes 
an arbitrarily chosen element of L^ oc (R 2 , M 2 ) such that J^i = Jj^ on L, the constant C being 
independent of this choice. 

If we choose an arbitrary Jo in A.\ and let the sum in (|7.32|) range over < i < N £ instead 
of 1 < i < N £ we obtain a measure P^ such that, by ()7.2ip . 



\\p {5) -p {6) \\<^ E 



1 JsTeJC„ ye 
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hence using Lemma 17. II we have d B (P {5) , P (6) ) < e and then d B (P,P {6) ) < Cefllog e| + 1). 

We now define the vector field j^ nt : M 2 — >■ M 2 by letting j^(x) = o mK 3^i,{x — x£) on 
Z = + L/ yJrriK-, for every .fT G /C n , < i < N £ and L G Then, for every L G we 

have J^ nt (x^ + •) = Om K Ji,L on L, therefore we may choose J^l = <r 1 / JTtx j™ t (x^ + •) in (|7,32p 
and then then we may summarize the above by writing 



(7.33) d B (P,P®) < Cefllog e| + 1), 



P (6) = T Ly 

\E'\ ^ 



K 



5 X K 



5fl, lint C?A. 



Note that since Jj l = outside L, we also have 



(7.34) 



•int 
J 11 



l<i<N £ 
LgCka 



0",, 



mp; J'i,L(,' 



curlj; nt = 2vr £ (<5 p - m K ) 



where A^- is a finite subset of the interior of K. The second equation is satisfied in the sense 
of distributions on IR 2 . 

Step 4 : Treating the boundary. Let E' := E' \ L)Keic„K . We let t G [0,^-v/n] denote arclength 
on dE' — where i is the length of dE — and s denote the distance to dE' , so that (t,s) 
is a C 1 coordinate system on {x G E' : d(x,(E') c ) < Cy/n}, if c > is small enough, 
since the boundary of E is C 1 by (|1.13|) . We let C% denote the curvilinear rectangle of 
points with coordinates in [0,t] x [0, CR e ], where R £ = q £ R £ and C is large enough so that 
E' C {x G E' : d(x,dE') < CR £ }, and define m(t) = Jc t nE' m °'' Since the distance of 
Uk£K.„K to a given x G dE' is between R £ and CqR £ from Lemma 17.51 and since mo is 
bounded above and below by ()1.15p . the derivative of t 1— > m{t) is between R £ /C and CR £ 
for some C > large enough. 
We let 



(7.35) 



Rr 



and choose = to, • • • , tk e = ^V™ to be such that 



m{ti) 



■rn 



We note that indeed tk £ = £\/n : Since the integral of mo' on each square K G K, n is an 
integer as well as the integral on E' , we have J g,, mo' G N and therefore m(£y/n) G N. 

Prom the above remark about the derivative of t — > m(t), we deduce that m< f^ 1 belongs 

to the interval [R £ /C, CR £ ] for some C > and then it is easy to deduce that if yjn is large 
enough compared to R £ then 



ni 



m{t 



z+i, 



m(U) G \R £ 2 /C,CR £ 2 



t l+1 -U G [R £ /C,CR £ ]. 



This means that the sidelengths of the curvilinear rectangle C ;+1 \C ; are comparable to R £ , 

2 

and that the number of points n\ to put there in is of order R £ . 

We may then include each of the sets K\ := E' n (Ct l+1 \ C^) in a ball B\ with radius 
in [R £ /C,CR £ ] and we may also choose a set of n/ points which are at distance at least 
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1/C from each other and the complement of K\. Let j\ — — V - *-/?, where H solves —AH = 
^^(SpeAi dp ~ m i) m &l an d d u H ■ t = on dBi, where 

ni 

mi = W\ lKr 

Then we have curlj; = 27r(^ pgA; 5 P — mi) in Bi and ji ■ r = on dBi and we claim that for 
any q > 1, 

(7.36) W(j h l Bl ) < C £ , UWl^kO < C £ ,g, 

where the constants do not depend on n, but do depend on e through R £ . This is proved 
by noting that these quantities are finite, and that a compactness argument shows that the 
bound is uniform for any choice of points which are at distance at least 1/C from each other 
and the complement of some K\ C B[, using for instance the explicit formulas for W in 
|LR| . Note that because the sets {K{\ and the rectangles {K} are disjoint, have measure 
between R e jC and CR e and diameter between R £ /C and CR E , we know that their overlap 
is bounded by a constant C independent of e, n. 

Step 5: Rectification of the weight. We rectify the weights mx, fxi\: For K £ fC n we let Hk 
solve —AHk = 27r(mo' — mx) on K and d u Hx = on dK. Similarly we let Hi solve the 
same equation with B[ replacing K and mi replacing mx ■ By elliptic regularity, we deduce for 
any q > 1 that \\V Hk\\li(k) (resp. \\V Hi\\ L q^ Bl ^) is bounded by C gi£ \\mo' — mK\\L°°(K) (resp. 
Cq )£ \\mo'— Since mo is C 1 we have |Vm,o'| < C/y'n, therefore ||mo'— fnK\\L°°(K) ^ 
CR £ /y/n, while ||mo' — mi\\ L c<>^ B ^ < C. We deduce that 

(7.37) \\VH K \\ Lq < ||V#,|| L , < C 9)6 . 
We let 

(7-38) jK=j™ t {K , 
and 

(7.39) j n = + ^ -V^tf* + ~ V±R i = E ^ " V±jf ^ + E " V± ^' 

KeK„ i=l KaK n 1=1 

A n = U KeK A K ^tl A «> 

where jx and -ff/c are set to outside K and similarly for ji, Hi outside Bi. Then curl j = 
^^(SpeAn ~ m o') m This completes the construction of j n . 

7.4 Estimating the energy 

Step 1: Energy estimate. We have 

W(j K ,l K )= J] W(a mK J itL (--x L ),L). 

0<i<N £ 
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From £LU) we find, letting L = y/m^{L - x L ), using ([7T29]) and \L\ = \K\/q 2 , that 
(7.40) W(j K ,l K )= Wia^J^l^J^lKlfe^Wia^J^ + CeY 



0<i<N £ 



We estimate the integral of | j n \ 2 on M 2 \ U p( z\ n B(p, 77). From (|7.39p . this integral involves 
on the one hand the square terms 

k 

(7.41) [ li/-V ± ^| 2 + E / Hk-^HkI 2 , 

where K, q = K \ U p£ \ n B(p,r]) and similarly for (Bi) v , and on the other hand the rectangle 
terms 

£ / tiK-V ± H K ).(j K ,-V ± H K/ )+ £ ■■■+ E ••• 

K,K'£K n Jk ^ K 'v l<l^i<k e KEK„ 

K^K' l<Z<fc £ 

We estimate the latter as follows: Since the rectangles in tC n do not overlap, the first sum is 
equal to zero. A nonzero rectangle term must involve some Bi, and moreover a given B\ can 
only be present in a number of terms bounded independently of n, £ because the overlap of 
the balls B\ and the rectangles K is bounded. Thus from (|7.35p we have at most C\fnjR e 
nonzero rectangle terms. Moreover, since the K^s are disjoint, and disjoint from the K's, 
in a rectangle term involving B\ n K the integral can be taken over K\K\, and in a term 
involving B\ n Bi it can be taken over [B\ fl £?i \ K^) U [B^ V\B\\ K\). 

In any case we use Holder's inequality and the bound — ^^Hi\\li{Bi\Ki) — Ce,<j f° r 
some q > 2, which follows from (|7.36p . f|7.37j) . together with the bound 

\\jl - V ^|Im'(B z )' ^ K ~ ^ ±Hk Wli'(K) - C e,qi 

which follows from (|7.40p . ()7.36p using Lemma 4.7 in [SSlj . to conclude that each rectangle 
term is bounded by C £ and then that their sum is 0(y/n), meaning a quantity bounded by a 
constant depending on e times y/n. 

The limit as rj — > of the terms in ()7.4ip is estimated as above by expanding the squares 
and using Holder's inequality with (|7.37|) . (|7.36|) . (|7.40p . together with the bound (|7.35p to 
show that 

Hm n^fe/ li/-V ± F / | 2 + V / \j K -V ± H K \ 2 +7T#A n log V ] 

< E W(j K ,lic) + 0(Vn). 

In view of the bound 0(\/n) for the rectangle terms and (|7.40p we find using (|7.29p that 

(7.42) W(j n ,l R 2)< E \K\^W(a mK Ji) + Cns + 0(^). 

KeK n Qe 

0<i<N £ 
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Step 2: We proceed to estimating W(j n , l^a). We have, using (|7.26p . (|7.25p . (|7.2ip . then the 
fact that m ' - m K < CR £ /^/n on K, then (T77TT|) with (OBI , then dTTTT]) and finally (17TTD . 
that 



V ~2~Wk (o~ m ic Ji) < \E'\ V P ( -= X ff* ) W(a mi< Ji) + \K\s 
i=l Qe i=l ^ n / 

< / WjfCo-mo^jJOciPCx^O + IA-l (4= +£ ) 

< V / (o-mo'^j) dP(a, j) + |# I ( -S= + Ce ) 

(7.43) = / WkU) dP(x,j) + \K\ (-^= + Ce) . 

Here we have used the fact that Wk is bounded below by some (negative) constant, a fact 
proved in [SSlJ that we use below several times. 

We proceed by estimating r^x- From (|7.26|> we deduce that 



VV™ J ~ \K\ \\E'\ \V™ 
and then it follows from (|7.29|) that 

N c „2IEV| 



£ _ 2| p/l / 7^ \ 

no* = <fe 2 - g < iV £ + x H £ c j . 

Summing over K 6 /C n , using the fact that 

(7.44) \E'\U Kn K\ <C £ y/n~ 
and then (|T2T]) . (fTTpl) . we find that 

^ ^^(a^ J ) < C|E'| x tf £ c ) + ±= +e^l < Cn + eY 

Summing (|7.43p with respect to K £ JC n and adding the above estimate we find, in view 
of (USD, flZaZP and (USD, that 

(7.45) W(j n , l Ra ) < n|E| / T^(j) dP(x, j) + Cn f e + V 

S^ep 5: Energy bound for {x\, . . . ,x n ). From (|7.45p . the constructed currents {j n } and points 
{A„,} n satisfy curl j n = 2ir \^2 p&An S p - m 'j in M 2 with #A n = n and 



(7.46) limsup rc^Jg?) < |£| / VKk-C?) dPfoi) + Cfe- 
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Now let {xi}i = {p/ \/n}p£A n be the points in A n in the initial scale, and let v n = <W 
Since A n C E' by construction, we have Supp(^ n ) C E and thus f ( dv = 0. Moreover, 
defining H' n by (|l,19p . we have that —AH' n = curl j n and we may thus write j n = —V^H n — 
Vf n for some function f n . But j n = outside of E, by construction, while H' n decays fast at 
infinity by its definition (|1 . 19[) and the fact that the right-hand side of (|1 . 19[) has integral 0. 
Letting U n = U pe A„B(p,7]), we first have 

f \jn\ 2 ~ f \jn + V/ n | 2 = -2 [ (j n + Vf n ) ■ V f n + f |V/ n | 2 . 

JB R \U V JB r \U v JB R \U n JB R \U V 

Since j n G L q loc for any q < 2 and since f n G W^(M. 2 ) for all q, the last two terms on 
the right-hand side converge as rj — > to the integrals over Br. Also integrating by parts, 
using the Jacobian structure and the decay of f n and j n , we have f B (jn + V/ n ) • V f n = 
j B —V^H' n ■ V f n — > as R — > +oo. Therefore, letting 77 — ^ then R — > +00 in the above 
yields 

W(j n , 1 R 2) - W(-V ± <, 1 R2 ) > J \Vf n \ 2 > 0. 

Together with (|7.46[) , we deduce in view of f)2. lj) that 

(7.47) limsup-(wn(xi,...,x n )-n 2 I(fio) + -logn) <— W K (j) dP(x,j) + Ce. 

n->oo U \ 2 J 7T J 



Step 4-' Existence of A n . We claim that if n is large enoug h and if j G Lf oc (IR 2 ,K 2 ) is such 
that 

(7.48) d p (j(yfcx + -),j^(V^x + •)) < m/2 
for any x G E\F for some set F satisfying \F\ < then 

(7.49) d B f£ 5 X ® 5 6ve j dx >j E 6 * ® 6 8^j^ dx) < Ce(|log e| + 1). 

This would follow immediately from Lemmas 17. 1| 17.21 and 17.31 if 0^ x j n nt belonged to some 
compact set independent of x ^ F and n. In our case we note that if x belongs to some 
L G Ck,i, where K G K. n and < i < N e , then 

j n nt (x + ■) = CJm K Ji,L(- +X- X L ). 

Moreover, since Jj G H e , from (|7.14|) it follows that if x — x^ ^ r(J,) / \jrax then / := 
o~m K Ji(- + x — x^) G G £ . If in addition, dist(x,dL) > rjoR £ , then we deduce from (|7.19|) that 
dp (j^ nt (x + •)>/) < J 7i/2 and d p (j(x + -),j') < ^l- Lemma I7T21 then yields ds(d J -int( a . + .), fy) < 
e and dB(Sjr x+ .-\,8j>) < e thus 

dB(5jtet( a . + .),(5 J -( a . + .)) < 2e. 

In view of Lemma 17.31 we find 



d B 



PI ^ ^ ® dx ' p| ^ . ® 5 e^j^ dx ) < Ce (i io g e i + !)» 
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where F is the union of F and of the union with respect to < i < N £ , K £ K, n and L £ vCj^i 
of ^. (x L + Y(Ji)/^mE), of ^{x £ L : dist(x,9L) < r? i? £ }, and of £ \ U^ n ^. It turns 

out that |F| < Ct]q if n is large enough, C being of course independent of e, and thus using 
Lemma ITdl we deduce (|7.49p . The claim is proved. 

To prove the existence of the set A n , we note that the currents Jj used in constructing 
j n nt depend on e but are independent on n. Then they are truncated to obtain J^x where 
the sidelengths of L are in [R £ /C,CR e ], i.e. in an interval independent of n. It follows at 
once that there exists 5 > such that the points in L may be perturbed by an amount 5 so 
that for every i, K and L £ Ck,% the perturbed jf L is at a distance at most r/i/4 of Ji t x, for 
every n. Then in view of (|7.44p and (|7,37p it follows that for n large enough the resulting jn 
will satisfy (|7.48p for x far enough from dE' , i.e. outside a set of proportion relative to \E'\ 
tending to as n — > oo. We deduce that jn satisfies (|7.49|) . hence if n is large enough 

d B {P 3 *,P) <Ce{\\og e| + l). 

The same reasoning implies that if we let {x\, . . . , x n } be the points in A n in original coordi- 
nates, then perturbing the points in A n by an amount 5 > small enough, i.e. perturbing the 
Xi's by an amount 5/y/n at most we obtain points yi such that w n {yi) < w n {xi) + e. Since 
the ordering of the points is irrelevant, we let S n denote the set of permutations of 1 ... n and 
define 

A n = {{yi, . . . ,y n ) :3a e S n \xi - y a ^\ < 5. 

Then, given rj > 0, from the previous discussion and choosing e > small enough we have 
for any n and any (yi, . . . ,y n ) £ A n that (|4,ip is satisfied and the existence of j n such that 
cml j n = 27r S y j — mo') and such that {Pj n } n satisfies (|4.2|) . 
This concludes the proof of Proposition 14.11 
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